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INTRODUCTION!  AND  SUMMARY 


A  military  need  exists  .  for  microwave  solid  state  device 
technology  with  complementary  analog  signal  processing  capabilities 
which  has  been  developed  over  the  past  ten  years  in  surface  acoustic 
waves  (SAW)  technology.  In  this  technology,  non-recurs i ve 
transversal  filters  have  been  realized  through  either  interdigital 
transducers  and/or  periodic  reflective  arrays  defined  on  the 
piezoelectric  substrates  in  which  the  surface  acoustic  waves  are 
propagated.  Utilizing  this  technology  simple  delay  lines,  complex 
matched  filters,  high  Q  resonators  and  chirp  filters,  have  been 
realized  in  the  VHP,  UHF ,  frequency  range.  Unfortunately,  extension 
of  this  technology  to  the  microwave  region  (f  >  IGHz)  has  been 
difficult.  At  these  frequencies,  the  surface  acoustic  wavelength  is 
less  than  lum  making  fabrication  difficult  and  the  propagation  loss 
becomes  excessive  (>100  db/ysec  at  10  GHz). 

Final  report  contract  AFOSR  80-8264  summarizes  the  results  of  a 
viable  technological  realization  of  this  goal  through  magnetostatic 
waves  (MSW)  propagating  at  microwave  frequencies  in  magnetically 
biased,  liquid  phase  epitaxial  films  of  Yttrium  Iron  Garnet 
(YIG)  grown  on  Gadolinium  Gallium  Garnet  (GGG) .  This  technology 
has  a  number  of  advantages;  low  loss  (<30  db/usec  at  10  GHz), 
tunable  by  bias  field  to  any  center  frequency  (1-40  GHz),  and 
finally  a  well  understood  and  characterized  wave  phenomena 
(MSW).  The  necessary  non-recursive  transducer  filtering  is 
achieved  through  a  mul t i -e 1 ement  transducer  and/or  the  interaction 
of  the  MSW  with  periodic  structures  in  or  on  the  YIG  film. 

This  report  summarizes  the  work  accomplished  toward  achieving 
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MSW  transversal  filtering.  First,  studies  of  a  variable  time  delay 
device  based  on  a  cascaded  MSSW  forward  wave  delay  line  and  an  MSBVW 
backward  wave  delay  line.  are  presented.  Next,  studies  of  MSSW 
(Magnetostatic  Surface  Wave)  transducers  composed  of  periodic  arrays 
of  narrow  shorted  microstrips  are  summarized.  A  3  port  model  for 
this  type  of  transducer  and  experimental  measurements  of  test 
devices  were  carried  out  to  verify  the  model.  A  synthesis  procedure 
has  been  developed  to  allow  synthesis  of  the  desired  response 
function.  Experimental  studies  of  synthesized  filters  has  resulted 
in  the  first  "true"  synthesized  MSW  transversal  filter.  Finally, 
studies  of  MSSW  and  MSFVW  delay  lines  at  Ku  Band  (26-40GHZ)  were 
carried  out  resulting  in  sample  delay  lines  with  less  than  16db 
insertion  loss  and  1  GHz  bandwidths. 

Pulse  compression  systems  utilizing  surface  acoustic  wave 
(SAW)  dispersive  delay  lines  have  been  under  intensive  investigation 
for  the  past  14  years.  These  studies  have  concentrated  primarily  on 
using  the  SAW  linear  dispersive  delay  line  for  downchirp  phase 
coding  with  quadratic  phase  error  less  that  one  degree  in  the 
expander  section  of  the  SAW  pulse  compression  system.  The  low  time- 
si  de  lobes  is  achieved  by  appropriate  amplitude  weighting  in  the 
compressor  section.  T ime-bandw i th  products  of  over  1000  with  greater 
than  40  dB  side  lobe  suppression  have  been  obtained  for  the  SAW  pulse 
compression  system.  The  typical  operating  center  frequency  for  a  SAW 
device  is  less  than  1  GHz  (due  to  the  physical  limit  of  technology). 
For  radar  systems  of  higher  operating  frequencies  (in  the  GHz 
region),  mixing  is  required  to  convert  the  radar  carrier  frequencies 
to  the  operating  frequencies  of  the  SAW  pulse  compression  subsystem. 
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A  new  technology  based  on  "slow"  magnetostatic  wave  (MSW) 
propagation  in  a  magnetically  biased  epitaxial  Yttrium  Iron  Garnet 
(YIG)  is  emerging  as  a  complementary  signal  processing  technology  to 
SAW  at  microwave  frequencies  (1-20  GHz).  Three  major  propagating 
modes  with  the  propagation  direction  in  the  film  plane  have  been 
used  in  device  applications.  These  three  modes  are  magnetostatic 
surface  wave  (MSSW)  mode  in  which  the  bias  field  is  perpendicular  to 
the  direction  of  the  wave  propagation  and  in  the  plane  of  the  film, 
magnetostatic  forward  volume  wave  (MSFVW)  mode  in  which  the  bias 
field  is  perpendicular  to  the  film  and  the  propagation  direction, 
and  magnetostat i c  backward  volume  wave  (MSBVW)  mode  in  which  the 
direction  of  the  bias  field  is  the  same  as  the  propagation 
direction.  These  modes  are  dispersive  and  characterized  by  a  limited 
propagation  passband  width  (.5  to  2.2  GHz)  with  the  center  frequency 
electronically  tunable  between  1  to  20  GHz  by  adjusting  the  bias 
field.  The  1-20  GHz  center  frequency  range  of  this  MSW  modes  makes 
signal  processing  possible  directly  at  radar  carrier  frequency  and 
the  typical  group  delays  of  the  modes  are  between  50  nsec/cm  to 
lOOOnsec/cm,  depending  on  the  YIG  film  thickness  and  the  bias 
magnetic  field  used.  This  dispersive  slow  wave  structure  of 
magnetostatic  wave  makes  the  investigation  of  MSW  devices  in  pulse 
compress  ion  logical. 
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2-  MAGMEnrOSTAT  I  C  WAVE 

D  I  SIXERS  I  ONI  CONITROL 

2.1  Linearly  Dispersive  Time-Delay  Control  of 
Magnetostatic  Surface  Wave  by  Variable  Ground  Plane 
Spac I ng 

2.1.1  I ntroduct I on 

The  propagation  character i st I cs  of  magnetostatic  waves  in 
epitaxial  Yttrium  Iron  Garnet  (YIG)  films  have  been  extensively 
investigated  both  experimentally  and  theoretically  for  a  number  of 
years  (1-4).  These  studies  have  shown  that  the  delay  versus 
freauency  characteristics  of  MSSW  delay  lines  of  a  YIG  film  spaced 
from  a  ground  plane  are  functions  of  the  saturation  magnetization, 
bias  field,  YIG  film  thickness,  and  ground  plane  spacing.  A  study 
of  how  different  parameters  affect  the  delay  characteristics  of 
surface  waves  indicated  that  dispersion  control  by  ground  plane 
spacing  is  most  effective  (4).  This  chapter  described  the 
thecret’cc'  and  experimenta’  results  on  the  theoretical  simulation 
and  experimental  evaluation  of  a  magnetostatic  surface  wave  variable 
ground  plane  spacing  delay  line  with  linear  group  delay. 

2.1.2  Theoretical  Simulation 

Figure  2.1.1  shows  the  variation  of  the  group  delay  with 
frequency  for  MSSW  delay  lines  with  ground  plane  spacing  (Ti)  and 
propagation  length  (Li)  as  adjusting  parameters.  The  YIG  film 
thickness  used  in  these  calculations  is  30  urn.  Examination  of  this 
figure  indicates  that  it  is  possible  to  achieve  a  linear  variation 
of  delay  by  adjusting  the  length  of  the  short  section  (Li)  of  the 
step  ground  plane  as  shown  In  Figure  2.1.2.  A  nine  step  ground  plane 
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structure  is  selected,  with  an  eight  step  ground  plane  in  the  middle 
of  this  MSSW  stepped  ground  plane  delay  line.  The  separation  of 
this  eight  step  ground  to  the  YIG  film  can  be  adjusted  by  a  pair  of 
set  screws.  The  ninth  step  in  this  MSSW  delay  is  fixed  by  the 
substrate  thickness  which  is  250  um.  The  propagation  length  can  be 
adjusted  by  photographically  changing  the  single  bar  transducer 
separation  between  it  and  this  edge  of  the  substrate.  The  length  of 
each  section  is  then  optimized  to  provide  a  minimum  deviation  from 
quadratic  phase  over  the  linearized  delay  bandwidth.  The  linear 
time  delay  characteristic  which  is  obtained  by  summing  the 
individual  time  delay  corresponding  to  each  section  of  the  MSSW  9 
stepped  ground  plane  structure,  is  shown  in  Figure  2.1.3.  A  linear 
regression  fit  is  also  shown  in  this  figure  which  indicates  a  I  inear 
delay  region  between  2.89  to  3.33  GHz  (Figure  2.1.4.).  The 
R.M.S.  deviation  from  quadratic  phase  over  this  440MHz  bandwidth  is 
8  degrees.  The  corresponding  time  bandwidth  product  (TBP)  is  near 
120  for  this  nine  stepped  ground  plane  MSSW  theoretical  simulated 
1  inear  delay  1  ine. 

2.1.3  Experimental  Evaluation 

Due  to  the  difficulty  and  long  waiting  period  for 
fabrication  of  the  eight  stepped  center  ground  plane,  a  continuously 
variable  ground  plane  is  used  to  provide  the  variation  of  ground 
plane  spacing  to  the  YIG  film  along  the  propagation  direction  of  the 
MSSW  delay  line.  The  separation  of  this  continuously  variable 
ground  plane  to  the  YIG  film  can  be  adjusted  by  changing  structure 
of  the  ground  plane  and  varying  the  distance  of  this  ground  to  the 
film  by  a  pair  of  set  screws. 
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Figure  2.1.1 


Variation  of  group  delay  with  frequency 
for  .MSSW  delay  lines  with  cround  plane 
spacing  (Ti)  and  propagation  distance 
(Li)  as  adjusting  parameters. 
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Ti( uM) 

Li ( cm ) 

300 

.18 

250 

.20 

225 

.15 

200 

.15 

175 

.23 

150 

.22 

135 

.20 

115 

.  10 

100 

.07 

Table  2.1.1  Ground  plane  spacings  and  lengths  for  a  nine- 

section  delay  line 
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TIME  DELAY  ( nanoseconds  ) 


FREQUENCY  (  G  Hz  ) 


Figure  2.1.3  The  summation  of  the  individual  time  delay  section 
in  the  9-step  ground  plane  structure,  with 
dimension  shown  in  table  2.1,  results  a  linear 
time  delay  characteristic  between  2.39  to  3.33 
GHz.  Thickness  of  YIG  is  30  uM. 


2.90  3.00  3.10  3.20  3.30  3.40 

FREQUENCY  (  GHz ) 

Figure  2.1.4  The  quadratic  phase  error  of  the  linear  delay 
region  of  the  9-3tep  ground  plane  structure 
The  R.M.S.  deviation  from  quadratic  phase  in  this 
bandwidth  is  8  degrees. 
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Figure  2.1.5  shows  a  sketch  of  the  MSSW  variable  ground  plane 
delay  line.  The  device  Is  In  a  "flipped"  configuration  with  the  YIG 
film  flipped  over  the  single  bar  mlcrostrlp  transducer.  A  pair  of 
shorted-c 1 rcu 1  ted  50  um-wide,  3  mm- long,  and  3  um-thlck  microstrip 
couplers  are  used  for  the  launching  and  receiving  of  the  microwave 
energy  for  MSSW  propagation  along  the  YIG  film.  A  25  um-thick  3mm- 
wide  YIG  film  and  15  um-thick  Gadolinium  Gallium  Garnet  (GGG) 
substrate  was  used  for  the  delay  measurement.  The  HP  8409B 
automatic  network  analyzer  was  used  to  measure  the  reflection  and 
transmission  parameters  of  the  device  under  test. 

The  group  delay  was  then  obtained  from  the  Sj2  phase  data  by 
calculation  of  phase  slope.  Since  this  Is  a  differentiation 
process,  the  resolution  is  determined  both  by  the  frequency  aperture 
(Af)  and  the  phase  detector  resolution  of  the  network  analyzer.  By 
increasing  the  frequency  aperture  and  the  corresponding  phase 
difference,  an  averaging  effect  is  performed  to  reduce  the  time 
delay  ripples.  A  30  MHz  frequency  aperture  Is  chosen  In  order  to 
optimize  the  calculated  time  delay  obtained  from  the  automatic 
network  analyzer  phase  measurement.  A  modification  of  the  control 
program  of  the  automatic  network  analyzer  Is  performed  to  compute 
the  absolute  phase  change  as  a  function  of  frequency  over  the 
frequency  bandwidth  of  Interest.  For  a  linear  delay  bandwidth,  the 
corresponding  absolute  phase  function  Is  a  quadratic  function  of 
frequency.  (For  a  constant  delay  bandwidth,  the  corresponding 
absolute  phase  function  Is  a  linear  function  of  frequency).  A  root 
mean  square  fit  is  then  performed  for  this  experimentally  determined 
phase  data.  The  deviation  from  quadratic  phase  defined  as  the 
experimental  phase  data  is  then  used  as  the  parameters  to  determine 
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A  15  STEP  3ROUNO  PLANE  MODEL  IS  USED  TO 
OBTAIN  THE  THEORETICAL  CASCADED  DELAY 
CHARACTERISTIC 
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Figure  2.1.5 


Sketch  of  continuously  variable 
delay  line 
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the  deviation  from  linear  delay  (quadratic  phase)  over  the  bandwidth 
of  interest  of  the  device  under  test. 

By  carefully  adjusting  the  separation  of  the  continuously 
variable  ground  plane  to  the  YIG  film,  linear  delay  of  bandwidth  up 
to  500  MHz  can  be  experimentally  obtained  with  phase  error  less  than 
20°.  The  slope  of  the  linear  region  can  also  be  adjusted  by 
changing  the  structure  of  the  continuously  variable  ground  plane. 
Figure  2.1.6  shows  the  experimental  delay  measurement  of  one  of  the 
devices  under  test.  A  linear  delay  of  500  MHz  bandwidth  (between 
2.62  to  3.12  GHz)  with  a  slope  of  15.3  nsec/ 100  MHz  was  observed. 
Figure  2.1.7  shows  the  quadratic  phase  error  of  this  device  in  the 
linear  delay  region.  The  corresponding  R.M.S.  phase  error  is  13°. 
The  insertion  loss  over  the  linear  region  is  between  1 5  to  20  dB 
(unmatched).  The  separation  of  this  continuously  variable  ground 
plane  to  the  YIG  is  roughly  determined  by  removing  the  "flipped" 
over  YIG  film  and  measuring  the  profile  of  the  MSSW  delay  line  along 
the  direction  of  propagation.  A  50um  air  gap  was  added  to  the 
measured  data  to  represent  the  air  gap  that  is  unavoidable  in  the 
flipped  configuration.  A  15-step  ground  plane  structure,  with 
dimensions  given  in  Figure  2.1.5,  was  then  used  to  approximate  the 
continuously  variable  ground  plane  structure  measured.  The  computed 
delay  characteristic  of  this  15-step  model  is  also  in  Figure  2.1.6. 
The  difference  between  this  computed  delay  characteristic  and  the 
experimental  linear  delay  characteristic  is  quite  large.  This  large 
discrepancy  is  due  to  the  error  associated  with  the  measurement  of 
the  bias  field  strength,  YIG  thickness,  and  separation  of  the  ground 
plane  to  the  YIG  film. 
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Figure  2.1.7  The  quadratic  phase  error  of  the  continuously 
variable  ground  plane  MSSW  delay  line 
The  R.M.S.  phase  error  is  13  degrees  over  the 
linear  delay  bandwidth. 


Applications  of  these  MSSW  linear  delay  devices  include  pulse 


compression  systems  (5,6),  compressive  receivers  and  variable  time 
delay  devices  utilizing  a  cascade  configuration  with  a  linearized 
backward  volume  wave  device  (7,8). 
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2.2  Electronically  Variable  Time  Delay  Using  Cascaded  Magnetostatic 
Wave  Delay  Lines 

2.2.1  I ntroduct i on 

Two  features  which  make  magnetostatic  wave  (MSW)  devices 
attractive  are  (a)  an  operating  frequency  range  above  1  GHz  and  (b) 
electronic  tunability.  Both  of  these  features  are  required  for 
electronically  variable  time  delay  components  used  in  phase  array 
radar  and  communication  systems.  The  basic  architecture  for  a 

variable  time  delay  device  based  on  MSW  technology  was  described 
originally  by  Sethares,  Owens  and  Smith  (7).  A  Schematic  of  their 
approach  is  illustrated  in  Figure  2.2.1.  Briefly,  two  MSW  delay 
lines,  one  biased  to  operate  In  the  backward  volume  wave  (MSBVW) 
mode,  and  the  other  biased  to  operate  in  the  surface  wave  (MSSW) 

mode,  are  cascaded  together.  The  dispersion  (l.e.,  the  delay  vs. 
frequency)  of  the  MSSW  device  has  a  positive  slope,  while  the  slope 
of  the  MSBVW  dispersion  is  negative.  If  the  absolute  magnitudes  of 
these  two  slopes  are  identical,  and  If  the  dispersions  are  linear, 
then  the  net  dispersion  of  the  cascaded  device  will  be  zero  over  a 

common  frequency  band.  If  the  frequency  domain  of  one  delay  line  is 

shifted  by  changing  the  strength  of  Its  applied  biasing  field,  then 
the  absolute  delay  of  the  cascaded  device  will  also  be  shifted. 
Thus,  different  values  of  constant  delay  can  be  obtained  by  simply 
changing  the  strength  of  the  bias  by  a  small  amount.  One  impediment 
to  the  realization  of  a  device  of  this  type  with  large  bandwidth  is 
that  in  general  the  dispersion  of  typical  MSW  devices  is  nonlinear. 
In  the  original  experiments  no  attempt  was  made  to  linearize  the 
individual  dispersions.  Even  so  the  results  were  generally 
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encouraging.  In  this  chapter  an  Innproved  version  of  the  original 
device  In  which  care  has  been  taken  to  linearize  the  dispersion  for 
both  the  surface  and  backward  modes  Is  reported. 

2.2.2  Theory 

From  the  study  of  the  dispersive  characteristics  for  MSSW  and 
MSBVW  delay  line.  It  Is  found  that  the  shape  of  the  dispersion  curve 
Is  strongly  dependent  on  the  thickness  of  the  film  and  the  spacing 
above  the  ground  plane.  It  has  been  demonstrated  (9)  that  the 
linearity  of  the  dispersion  can  be  substantially  Improved  by 
carefully  adjusting  the  spacing  between  the  ground  plane  and  the  YIG 
film.  Thus,  the  ground  spacing  technique  was  Incorporated  Into  the 
design  of  the  delay  line  components  of  the  cascaded  device. 

A  linear  MSBVW  dispersion  curve  can  be  obtained  In  a  straight 
forward  manner  by  using  a  35  urn  thick  layer  of  YIG  spaced  250  um 
above  the  ground  plane.  This  spacing  Is  very  close  to  the  standard 
for  MSW  delay  1 Ines,  and  no  special  fabrication  techniques  are 
required,  other  than  some  polishing  of  the  dielectric  spacer.  It 
requires  a  much  closer  ground  plane  spacing  to  achieve  the  same 
results  with  the  MSSW  mode  and  past  experience  has  shown  that  ground 
current  losses  Increase  as  the  spacing  decreases. 

To  circumvent  this  problem  a  variable  ground  plane  structure 
of  the  type  reported  by  Chang,  Owens,  and  Carter  (10)  was  selected 
for  the  MSSW  device.  The  essence  of  this  approach  Is  that  the  YIG 
delay  line  Is  held  above  a  ground  plane  whose  distance  from  the  YIG 
film  Is  a  function  of  the  path  length.  A  schematic  of  the  technique 
Is  shown  In  Figure  2.2.1,  with  experimental  results  shown  In  Figure 
2.2.2  and  Table  2.2.1.  The  advantage  of  this  approach  Is  that  In 
principle  optimum  dispersion  control  can  be  achieved  and  the  minimum 
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Table  2.2.1  The  effect  on  changing  the  slopes  of  the  MSBVW  linear 
delay  line  by  adjusting  the  bias  field 


21 


spacing  required  is  greater  than  in  the  uniform  spacing  case. 

For  both  modes  of  propagation  the  figure  of  merit  is  defined 
in  terms  of  the  phase  deviation  from  the  appropriate  theoretical 
phase  vs.  frequency  curve.  If  the  delay  time  is  a  linear  function 
of  the  frequency  then  the  phase  dependence  on  frequency  is 
quadratic.  On  the  other  hand,  if  the  delay  time  is  constant  over 
the  frequency  band,  then  the  phase  is  a  linear  function  of 
frequency.  In  principle,  it  is  possible  to  achieve  close  to  zero 
deviation  from  quadratic  phase  over  a  bandwidth  of  500  MHz  centered 
at  3  GHz  using  the  ground  plane  spacing  techniques  described  in 
Chapter  2 . 1 

2.2,3  Experimental  Results 

A  cascaded  device  was  prepared  using  a  MSSW  delay  line  with  a 
variable  ground  plane  and  a  MS8VW  delay  line  with  a  YIG  film 
thickness  of  35  um  and  ground  plane  spacing  of  254  um.  The 
ch  racteristic  of  the  MSBVW  delay  line  is  shown  in  Figure  2.2.2. 
The  variable  ground  plane  was  prepared  by  polishing  a  ground  plane 
into  a  smooth  curve.  The  ground  plane  curve  was  optimized  by 
alternating  vs.  frequency  measurements  with  block  sanding  until  the 
desired  dispersion  curve  was  obtained. 

The  delay  line  material  was  prepared  by  growing  YIG  films  on 
the  I  inch  diameter  Gadolinium  Gallium  Garnet  (GGG)  wafers  using  the 
Liquic  Phase  Epitaxy  (LPE)  method.  The  YIG  thickness  for  the  MSSW 
delay  line  was  25  um,  while  the  thickness  of  the  MSBVW  delay  line 
YIG  was  35  um.  The  wafers  were  cut  into  bars  for  delay  line  use, 
and  the  ends  of  these  bars  were  beveled  to  reduce  reflections. 

The  individual  delay  lines  were  experimentally  evaluated  with 
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an  automatic  network  analyzer.  Large  laboratory  electromagnets  were 
used  to  supply  the  necessary  d.c.  bias.  In  the  final  configuration 
both  the  MSSW  delay  line  and  the  MSBVW  delay  line  exhibited  R.M.S. 
phase  deviation  from  quadratic  of  about  13°  over  approximately 
500  MHz  bandwidth  centered  at  3GHz.  Experimental  plots  from  the 
network  analyzer  of  delay  time  vs.  frequency  for  both  MSSW  and  MSBVW 
devices  are  given  in  Figure  2.2.3.  The  slopes  of  these  curves  are 
9.6  nsec/100  MHz  and  -9.7  nsec/100  MHz,  respectively. 

It  is  found  that  by  adjusting  the  d.c»  bias  field  of  the 
variable  ground  plane  MSSW  delay  line,  the  dispersion  characteristic 
of  the  linear  delay  bandwidth  becomes  nonlinear.  In  fact,  the 
deviation  from  quadratic  phase  over  the  frequency  band  of  interest 
increases  substantially  by  changing  the  bias  field.  For  the 
backward  volume  wave  delay  line,  the  linear  bandwidth  can  be  changed 
by  adjusting  the  bias  field.  The  dispersion  characteristic  over  the 
bandwidth  of  interest  remains  linear  with  a  slight  adjustment  in  the 
slope.  Table  2.2.1  gives  the  effect  of  changing  the  d.c.  bias  of 
the  MSBVW  delay  line  on  the  slopes  of  the  linearized  frequency  band. 
Because  of  these  experimental  results,  the  bias  field  for  the  MSSW 
linear  delay  line  is  fixed  at  the  same  value,  whereas  the  bias  field 
of  the  MSBVW  linear  delay  line  can  be  changed  to  provide  adjustment 
for  the  time  delay  in  the  cascaded  device.  These  two  MSW  delay 
lines  are  then  cascaded  together  to  provide  the  variable  time  delay 
exper iments . 

Delay  vs.  frequency  curves  experimentally  obtained  from  the 
cascaded  MSSW  and  MSBVW  device  are  shown  in  Figure  2.2.4.  To  obtain 
these  curves  the  bias  on  the  MSSW  delay  line  was  held  constant  at 
410  Gauss  while  the  bias  applied  to  the  MSBVW  delay  line  was  varied 
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DELAY  (  pSEC) 


Figure  2.2.3  MSSW  and  MSBVW  dispersion  curve  and  their  sum 
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Figure  2.2.4  Dispersion  curves  for  cascaded  delay  line  as  a 
function  of  bias  field 
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from  540  to  650  Gauss.  The  data  obtained  are  summer i zed  in  Table 
2.2.2.  The  maximum  delay  differential  obtained  was  42nsec  with  an 
applied  bias  field  of  110  Gauss.  The  bandwidth  of  this  variable 
delay  line  is  between  2.9  and  3.15  GHz,  while  the  R.M.S.  phase 
deviation  from  linearity  ranged  from  6.4®  to  12.9®.  These  values 
are  quite  acceptable  for  a  number  of  phased  array  radar/communi ca¬ 
tions  applications  and  represent  a  significant  advancement  in  the 
state-of-the-art  for  MSW  devices. 

There  is  one  problem  with  the  cascaded  device  as  it  stands 
now.  The  insertion  loss  is  high,  approximately  35  dB  across  the 
band.  This  high  value  is  a  result  of  the  device  not  being  tuned.  The 
loss  could  be  reduced  considerably  with  proper  matching  of  the 
transducers . 

2.2.4  Summary 

An  electronically  variable  time  delay  device  comprised  of  two 
cascaded  MSW  delay  lines  has  been  fabricated  and  evaluated,  and  the 
operating  characteristics  are  substantially  better  than  previously 
reported  embodiments  of  this  device.  The  improvements  are  due  mainly 
to  the  1 inearizat ion  of  the  two  component  dispersions.  With  a 
reduction  in  insertion  loss  and  compact  packaging,  a  device  of  this 
type  will  be  ready  for  service  in  phased  array  systems. 
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Experimental  results  -For  the  MSW  Variable 

Time  de 1  ay  lines 


T  C  nsec ) 


R . M . S .  (degrees : 


540.0 

149.6 

9.75 

550.0 

154.7 

12.90 

562.5 

•153.7 

•11.10 

575.0 

164.7 

10.60 

587 . 5 

163.3 

9 . 50 

600.0 

173.8 

9.10 

6  12.5 

178.0 

6.40 

625.0 

102.4 

6.00 

637.5 

106.7 

6.40 

650.0 

19  1.4 

12.20 

Bandwidth!  2  .S 

TO  3.15  GHz 

Table  2.2.2  Experimental  results  for  the  MSW  Variable  Time 
delay  lines 
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2.3  Passive  Magnetostatic  Wave  Pulse  Compression  Loop 
2.3.1  I ntroduct i on 

A  MSW  pulse  compression  loop  using  active  expended  chirp 
signal  generation  and  magnetostaic  surface  wave  (MSSW)  delay  line 
with  a  nonlinear  delay  characteristic  as  the  compression  filter,  was 
first  demonstrated  by  K.W.  Reed,  et.al  of  this  laboratory  [11],  The 
"down  chirp"  signal  was  generated  by  a  voltage  controlled  oscillator 
driven  by  a  programmable  sweep  voltage  with  a  matched  nonlinear 
phase  characteristic  appropriate  to  the  intrinsic  non  linearity  of 
MSSW  delay  line.  Phase  predistortion  techniques  were  used  to  achieve 
a  24  dB  side  lobe  ratio  with  a  time-bandwidth  product  of  30. 

In  this  study  a  matched  pair  of  quadratic  phase  MSW  delay 
lines  with  a  center  frequency  of  3.2  GHz  and  3  dB  bandwidth  of  320 
MHz  have  been  used  as  expander  and  compressor  devices  for  pulse 
compression  loop.  A  schematic  diagram  of  the  MSW  pulse  compression 
loop  is  shown  in  Figure  2.3.1.  Appropriate  amplitude  weighting  by 
using  loop  transducers  in  the  MSSW  device  have  provided  30  dB 
side  lobe  suppression. 
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BLOCK  OIAGKMf  Of  MSH  PULSE  CGNPflCSSlON  LOOP 


Figure  2.3.1  Block  diagram  of  MSW  pulse  Compression  loop 
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2.3.2  Passive  Generation  Of  The  Down  Chirp  Coded  Signal 

A  simple  linear  backward  volume  wave  delay  line,  driver,  bv  a 
fast  rise  pulse  generation,  was  used  for  the  passive  generat l on  of 
the  FM  down  chirp  signal.  The  linear  MSBVW  dispersion  character i st i c 
was  obtained  in  a  straight  forward  manner  by  us  i  rig  a  50  urn  thicK 
layer  of  YIG  film  spaced  21'iO  um  above  the  ground  plane.  Tiie  1  i  near- 
dispersive  curve  of  this  delay  line  is  showri  in  Figure  Z.2.3.  The 
slope  of  this  linear  delay  line  is  -9.  i  nsec/lOO  MHz.  w'tl"'  deviation 
from  quadratic  phase  iess  than  10  degrees  over  500  MHz  bandwidth  at 
3  GHz  center  frequency.  The  insertion  loss  of  this  delay  1 ine  Is  Il¬ 
ls  dB  over  the  500  MHz  b^jndwidth. 

2.3.3  MSSW  Linear  Matching  Filter 

The  "matched"  filter  required  in  this  pulse  compression  i ooo 
is  a  de 1  ay  line  with  a  t i tne- i nverted  replica  of  the  expanded  signal 
as  its  impulse  response.  Since  a  linear  MSBVW  delav  witn  a  "down- 
cfi  i  rp"  quadr  at  i  c  phase  i.haracter  i  st  i  c  i  s  used  to  generate  tne  down 
chirp  signal  ,  a  I  inear  delay  with  opposite  constant  delav  ;  i  jpe  O'ver 
Che  same  frequency  band  must  be  empioved  as  the  matched  fi 'ter.  The 
vai  iable  ground  plane  spacing  techniuue  described.  i  r:  Chapter  1.  1. 
has  been  used  to  achieve  the  required  quadratic  cr  asa?  chir--" 

char  acter  i  ot  i  cs  in  a  MSSW  delav  i  ine.  I.  1  0  ]  .  'he  siopo  -  c  cns 
dispersion  charact er i ^ t i c ,  shown  in  Figure  1.1.3,  is  9.7  nsec, 10c 
MHz  with  a  R.M.S.  phase  ei  ror  less  tnan  11  deqre-s  o-er  .5  -.00  .'■'nz 

bandwidth  at  3  GHz  center  freuuency. 

2.3.4  Experimental  Results 

Initial  experimental  results,  using  s i ng , e  Par  trancducji  s  in 
tfie  MSSW  delay  !  ine,  sirowed  on  i  v  10  JB  o*^  s  :  de  I  jbe  sunor  es  i  i  on  . 
Th.is  is  a  result  of  incorrect  ampl  ;  tude  weighting  arn.J  prase 
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distortion  at  the  band  edge.  The  3  dB  pulse  width  is  approximately 
5.5  nsec  with  a  time  bandwidth  product  of  24.  An  improved  MSSW 
device  was  designed  with  a  loop  transducer  to  achieve  a  better 
weighting  function  characteristic  for  the  MSSW  pulse  response 
envelope.  The  deviation  from  the  Hamming  weighting  function  was  +  or 
-2.5  dB  over  the  320  MHz  bandwidth  at  3.2  GHz  center  frequency  for 
a  loop  transducer  of  250  urn  loop  spacing.  The  quadratic  phase 
behavior  was  matched  by  adjusting  the  ground  plane  according  to 
chapter  2.  A  30  dB  peak  side  lobe  ratio  was  obtained  by  using  this 
loop  transducer  MSSW  linear  delay  as  the  compressor  in  the  MSW  pulse 
compression  loop  with  no  phase  predistortion.  The  characteristic  of 
the  MSBVW  linear  delay  line  down  chirp  filter,  used  as  the  impluse 
expander  in  the  MSW  compression  loop,  is  summarized  in  Table  2.3.1. 
The  characteristic  of  the  MSSW  linear  delay  line,  used  as  the 
compressor  in  the  MSW  compression  loop,  is  summarized  in  Table 
2.3.2.  The  recompressed  signal  obtained  by  this  pair  of  MSW  matching 
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THE  CHARACTERISTIC  OF  MSBVW  LINEAR 
DELAY  LINE  DOWN  CHIRP  FILTER 

YIG  THICKNESS;  50  um 
SPACING  TO  GROUND:  250  um 
BANDWIDTH: 500  MHZ 
CENTER  FREQUENCY:  3.2  GHz 
BIAS  FIELD:  590  Gauss 
TRANSDUCER  SEPARATION:  1  cm 
INSERTION  LOSS:  12  TO  15  dB 
DEVIATION  FROM  QUADRATIC  PHASE:  10° 

SLOPE  OF  LINEAR  DELAY;  -9.1  nsec/100  MHz 


Table  2.3.1  The  characteristic  of  the  MSBVW  linear  down 

expansion  filter 


ch  I  rp 
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THE  CHARACTERISTIC  OF  THE  MSSW  LINEAR 
DELAY  LINE  COMPRESSION  FILTER 

VARIABLE  GROUND  PLANE  SPACING 

YIG  THICKNESS:  25  um 

BANDWIDTH:  320  MHz 

CENTER  FREQUENCY:  3.2  GHz 

BIAS  FIELD:  410  Gauss 

DEVIATION  FROM  QUADRATIC  PHASE:  12° 

SLOPE  OF  THE  LINEAR  DELAY:  9.7  nsec/ 100  MHz 
WEIGHTING:  loop  transducer 

250  uM  loop  spacing 


Table  2.3.2  The  characteristic  of  the  MSSW  linear  compression  filter 
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filters  is  shown  in  Figure  2.3.2  with  8  nsec  3  dB  pulse  width  and  a 
time  bandwidth  product  of  12. 


2.3.5  Cone  1  us i on 

A  pulse  compression  loop  using  a  pair  of  matched  MSW  delay 
lines  with  passive  (impulse)  generation  has  been  demonstrated  with 
-30  dB  time  side  lobes  and  a  time  bandwidth.  Finally,  a  comparison 
of  the  SAW  and  MSW  pulse  compression  loop  is  summarized  in  Table 
2.3.3.  As  shown  in  this  table  the  compression  gain  of  MSW 
compression  loop  is  modest  by  SAW  standards.  However,  significant 
potential  for  increase  in  gain  exists,  as  well  as  good  potential  for 
better  side  lobe  suppression. 
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Figure  2.3.2  Experimentally  recompressed  pulse  from  loop 
transducer  device 
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SAW 


MSW 


Tab  I  e 


IF  freq  (<1  GHz) 

-3  dB  pulse  width 
8  to  1000  nsec 

we i ght  f  ng 

using  interdigital 

transducer 

TBP  50  to  1500 


side  lobe  levels 
30  to  45  dB 


RF  freq  (1-20  GHz) 

-3  dB  pulse  width 
2  to  10  nsec 

we i ght i ng 

approximate  by  loop 
transducer 

TBP  12  to  30 
(experimental ) 

side  lobe  levels 
24  to  30  dB 


2.3.3  A  comparison  of  SAW  and  MSW  pulse  compression  loop 
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2.-4  Magne-tos-t^-tn  f  c  Si_iT-f=-ace  Wav's 
Ma-tcrI-i  1  rig  IMa‘tiwor“l< 

Impedance  matching  provides  a  method  of  maximum  power 
transfer  between  the  signal  source  and  the  load.  Usually,  a 
quarter  wavelength  microstrfp  line-section  Is  required  for 
Impedance  matching  when  the  value  of  the  source  resistance  is  quite 
different  from  that  of  the  load  Impedance.  However,  the  quarter 
wavelength  is  bandwidth  limited.  From  early  work  of  microstrip 
transducer  theory  and  experiments.  It  was  known  that  the  Input 
Impedance  of  a  single  bar  transducer,  used  as  a  microwave  energy 
coupler  In  the  "flipped”  configuration  with  the  YIG  film  flipped 
over  the  single  bar  microstrip  transducer,  depends  on  the  length  and 
width  of  the  single  bar  mIcrostrIp  and  the  thickness  of  the  YIG 
film  [12].  Experimentally,  the  coupling  between  electromagnetic 
waves  and  MSW  Is  strong,  and  the  resultant  Input  Impedance  of  a 
single  bar  MSW  mIcrostrIp  transducer  of  length  in  the  order  of 
•millimeters  is  In  the  range  of  50  Ohms.  Theoretically,  by  adjusting 
the  length  of  the  single  bar  transducer,  the  real  (resistance)  part 
of  the  input  Impedance  can  be  selected,  so  that  it  is  close  to  the 
source  resistance  of  50  Ohms.  Therefore,  by  adjusting  the  length  of 
the  single  bar  transducer.  Input  Impedance  matching  can  be 
performed.  The  length  of  the  mIcrostrIp  Is  then  adjusted 
accordingly  to  obtain,  to  first  order,  the  Input  matching  circuit. 


Usua 1 1 y , 

the 

Imaginary 

( reactance) 

part  of 

thi  3 

ml  crostr 1 p 

transducer 

1  s 

1 nduct 1 ve ; 

one  way  to 

compensate 

th  1  s 

1 nduct 1 ve 

reactance 

1  s  to 

Introduce  a 

series  Input 

capac 1  tor . 

Another  way  1  s 

to  use  one  short  section  of  mIcrostrIp  transmission  line  to  rotate 
the  Input  impedance  to  the  real  axis  of  the  Smith  chart  and  a 
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quarter  wavelength  line-section  of  appropriate  characteristic 
impedance  to  match  this  inductive  reactance  to  the  source  impedance 
that  i s  usua 1  I y  50  Ohms . 

The  objective  of  this  study  is  to  reduce  the  insertion  loss 
associated  with  the  mismatch  in  the  continuously  variable  ground 
plane  magnetostatic  surface  wave  linear  delay  device.  This  device 
was  50  um-wide  by  3  mm- long,  had  a  3  um  thick  microstrip 
transducer  and  had  a  25  ym  YIG  flipped  over  on  it.  It  was 
measured  with  the  external  bias  field  of  410  Gauss  with  the  output 
short  circuited  in  the  MSSW  mode.  The  HP  8409B  automatic  network 
analyzer  was  used  to  measure  the  reflection  parameters  of  the 
device.  The  reference  plane  chosen  in  this  measurement  was  in  the 
end  of  the  50  Ohms  cable  with  a  SMA  male  connector.  The  rotation 
of  the  input  impedance  data  due  to  the  electronic  length  of  the  SMA 
connector  and  the  50  Ohm  microstrip  line  is  accomodated  by  using  a 
NEG  statement  in  the  COMPACT  analysis  program.  The  corrected  input 
impedance  is  shown  in  Figure  2.4.1.  The  bandwidth  of  interest  in 
this  case  is  from  2.9  to  3.15  GHz  or  only  250  MHz.  The  resistance 
of  this  set  of  data  from  2.9  to  3.15  GHz  is  close  to  50  Ohms.  One 
can  rT>atch  the  reactance  of  this  set  of  data  to  50  Ohms  by  a  series 
capacitor  of  1.9  pF  or  by  using  a  50  Ohm  30  degree  microstrip  line 
to  rotate  the  parameters  to  the  resistance  axis  on  the  Smith  Chart 
and  a  65  Ohm  quarter  wavelength  transformer  to  match  the 
parameters  to  50  Ohms.  Figure  2.4.2  shows  the  theoretical 
predicted  reflection  coefficients  by  COMPACT  by  using  the  50  Ohm 
30  degree  microstrip  line  and  the  65  Ohm  90  degree  microstrip  as 
the  matching  network.  An  experimental  matching  network  using  the 
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Figure  2.i.l  The  reflection  coefficients  of  the  MSSW  devic 
a  shorted-circuited  50  um-wide,  3  nun- long,  an 
thick  raicrostrip  transducer  in  the  "flipped" 
configuration 
Thickness  of  YIG  is  25  ym 
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Figure  2.4,2  Theoretical  predicted  results  by  COMPACT  on  the 

ref  lectio.”,  coefficients  by  using  a  50  Ohm  30  degrees 
microstrip  line  and  a  65  Oh.m  90  degree  microstrip 
line  as  the  matching  network 


40 


50  Ohm  30  degree  microstrlp  Mne  and  the  65  Ohm  90  degree 
microstrip  transformer  was  built  and  tested.  Figure  2.4.3  shows  the 
measured  data  of  this  device.  This  is  very  close  to  the  theoretical 
predicted  data  shown  in  Fig  2.4.2. 

Experimentally,  one  finds  that  the  reflection  coefficients 
are  the  same  for  devices  with  or  without  the  continuously  variable 
plane  structure,  but  the  insertion  loss  for  the  device  with  a 
continuously  variable  ground  plane  is  higher  due  to  ground  current 
loss.  Therefore,  the  matching  network  should  work  for  a  device  with 
or  without  a  ground  plane.  By  matching  the  input  impedance,  the 
loss  due  to  reflections  is  minimized.  Figure  2.4.4  shows  the 
insertion  loss  of  the  MSSW  delay  line  without  the  variable  ground 
plane  structure  and  without  the  matching  network  from  2.8  to  3.3 
GHz;  the  maximum  insertion  loss  from  2.9  to  3.15  GHz  is  about  15  dB 
whereas  the  minimum  insertion  loss  in  this  frequency  band  is  11  dB. 
Figure  2.4.5  shows  the  improvement  due  to  the  matching  network  for 
the  MSSW  delay  line  without  the  center  variable  plane  structure;  the 
maximum  insertion  loss  is  12  dB  and  the  minimum  value  is  8  dB.  The 
overall  improvement  is  about  3  dB  across  the  bandwidth  of  interest. 
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Figure  2.4.3  The  experimental  measured  data  on  the  reflection 
coefficients  of  the  matching  network  by  using  a 
50  ohms  30  degrees  microstrip  line  and  a  65  ohms 
ouarter  wavelength  microstrip  transformer 
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Figure  2.4.4  The  insertion  loss  of  the  non-matched  MSSW  delay 
line  2.8  to  3.3  GHz 


Figure  2.4.5  The  insertion  loss  of  the  matched  MSSW  delay  line 
from  2.9  to  3.15  GHz 
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3.  MSW  TRAIMSVERSAL  RIL-TER 

3.1  MSSW  Transversal  Filter  Based  on  Open  Circuit  Transducer 
Arrays 

Although  some  limited  types  of  bandpass  shaping  are  possible 
using  multi  bar  transducers  or  reflective  arrays,  the  performance  of 
these  devices  exhibit  high  insertion  loss.  The  complexity  of  the 
theoretical  multi  bar  transducer  response  using  the  reflective  array 
techniques  (etched  metal  film  patterns,  ion  implanted  zone,  etc.)  is 
the  main  reason  for  this  group  to  seek  alternative  techniques. 

The  first  approach  was  to  investigate  the  practicality  of  the 
design  of  transversal  filters  by  current  weighting  of  multibar 
transducers.  To  do  so,  two  sets  of  single  bar,  3mm  I/O  transducers 
were  built.  One  set  was  made  using  the  standard  photo  I i thography 
and  gold  up-plating  technique  which  resulted  in  a  circuit  pattern 
with  gold  thickness  of  about  3  um.  The  ohmic  resistance  of  the 
transducer  was  approximately  0.5  Ohms.  For  the  other  set,  the 
transducers  were  covered  with  photoresist  before  gold  plating  the 
circuit  to  prevent  it  from  plating,  so  only  the  50  Ohm  transmission 
line  was  gold  plated  to  a  thickness  of  3  ym.  The  transducers  had  a 
thickness  of  about  800  Angstroms  and  the  ohmic  resistance  was  about 
200  Ohms.  The  theoretical  and  experimental  insertion  loss  vs. 
frequency  response  of  the  two  devices  are  shown  in  Figure  3.1.  The 
same  procedure  was  repeated  for  a  two  parallel  I/O  bar  and  a  four 
parallel  I/O  bar  delay  lines  with  transducer  spacing  of  300  um  for 
all  of  them.  The  results  are  shown  in  Figures  3.2  and  3.3. 
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Figure  3. 1  The  theoretical  frequency  response  of  a  singlebar  I/O  (TOP)  and 
the  experimental  result  for  3  micron  thick  transducer  (MIDDLE) 
and  800  angstrom  thick  transducer  (BOTTOM). 
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Figure  3.2  The  theoretical  firequency  response  of  a  2  parallel  Ibars  I/O  delay  line 
(TOP)  and  the  experimental  results  for  3  micro  thick  transducers 
(MIDDLE)  and  800  angstrom  thick  transducers  (BOTTOM). 
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Figure  3.3  The  theoretical  frequency  response  of  a  4  parallel  bars  VO  delay  line 

(TOP)  and  experimental  results  for  3  micron  thick  transducer  (MIDDLE) 
and  800  angstrom  thick  transducer  (BOTTOM). 


The  next  step  was  to  see  If  some  kind  of  current  weighting 

could  be  done  by  having  thin  gold  In  a  portion  of  the  transducer 

while  the  rest  of  the  transducer  length  had  thick  gold.  As  an 
example,  consider  two  delay  lines,  each  having  a  single  bar  in  the 

input  and  the  output  with  one  device  having  thick  transducers  and 

the  other  with  thin  transducers.  If  the  average  current  In  the 

thick  transducer  Is  I  and  the  average  current  in  the  thin  transducer 

Is  (where  obviously  I  >  Iq)»  then  the  minimum  current  weighting 

factor  is  equal  to  I  /I=a  where  0<a< 1 .  If  the  thickness  in  the 

o 

thin  transducer  is  very  small,  then  0«a<l.  Therefore  any  current 
weighting  could  be  obtained  by  dividing  the  transducer  into  two 
specified  regions  of  thick  and  thin  films.  Three  sets  of  5  bar 
I/O  transducers  were  fabricated.  One  set  had  thick  gold  on  the 
entire  length  of  each  transducer,  and  was  used  as  the  control 
sample.  The  responses  of  the  other  devices  were  compared  with 
this  control.  In  the  other  two  sets,  a  portion  of  each  transducer 
was  made  to  have  a  thickness  of  about  800  Angstroms  (shown  as  the 
white  area  in  figure  3.4),  and  the  remaining  length  had  a  thickness 
of  about  3  urn  (black  region). 

It  should  be  pointed  out  that  at  this  point  no  attempt  was  made 
to  precisely  control  the  ratio  of  the  thick  and  thin  gold  regions. 
The  main  interest  was  to  experimentally  Investigate  the  effect  of 
current  weighting  (although  unknown),  on  the  frequency  response  of 
the  multibar  transducers.  The  Insertion  loss  responses  of  the  three 
devices  are  shown  in  Figure  3.4.  As  can  be  seen  a  dramatic  change 
can  occur  if  the  current  Is  not  distributed  evenly  in  each  finger  of 
a  multibar  transducer  delay  line.  These  results  indicate  that  the 
frequency  response  can  be  changed  with  proper  current  weighting 
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technique. 

The  theoretical  approach  for  designing  a  filter  Is  given  In 
Appendix  1.  The  set  of  computer  programs  used  for  synthesizing  a 
filter  Is  printed  In  Appendix  2.  All  the  programs  were  written  In 
FORTRAN  and  were  run  on  a  DEC-20  system.  The  first  program  gives 
the  linear  Insertion  loss  vs.  frequency  response  (total  of  152 
points)  for  a  single  bar  delay  line  with  zero  path  length  in  the 
MSSW  mode.  The  second  program,  [13]  calculates  the  complex  K  values 
for  the  frequency  range  determined  from  running  the  first  program. 
Finally,  the  synthesis  program  was  run  which  used  linear 
frequency  vs.  Insertion  loss  data  and  complex  K  data  from  the  two 
previous  programs  to  generate  normalized  current  weighing  factors 
for  an  N  bar  transducer.  The  shape  of  the  filter  was  In  this  case 
predetermined  to  have  a  square  passband.  In  this  program,  one  can 
change  the  shape  of  the  passband,  the  bandwidth,  the  center 
frequency,  the  side  lobe  suppression  and  the  number  of  transducers. 
As  an  example,  the  program  was  run  for  synthesizing  a  bandpass 
filter  centered  at  2.8  GHz.  This  showed  the  effect  of  the  passband 
and  also  the  number  of  transducers  on  the  frequency  response  of  the 
filter.  The  passbands  being  considered  are  100  MHz,  150  MHz  and  200 
MHz.  For  each  case  the  number  of  transducers  varied  as  5,  7,  9,  11, 
13,  15  and  25.  From  the  results  shown  In  Figure  3.5,  It  Is  clear 
that  for  a  given  bandwidth  there  Is  an  optimum  number  of 
transducers.  Increasing  the  number  of  bars  beyond  this  point  will 
not  Improve  the  response  of  the  filter  any  further.  Figure  3.6  Is 
the  theoretical  response  of  a  2 1  bar  transducer  programmed  to  be  a 
100  MHz  bandpass  filter  centered  at  2.8  GHz  with  sidelobe 
suppression  of  100  dB.  As  can  be  seen  the  maximum  theoretical 
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er  of  transducers  are  varied.  These  sets  of  plots  Indicate 
re  is  an  optimum  number  of  the  bars  for  a  given  bandwidth. 


NUMBER  OF  BAR 


requency  response  and  its  expanded  form  of 
tslgned  CO  have  a  bandwidth  of  100  MHz, 

Iz  and  CO  have  a  sidelobe  suppression  of  100 


side  lobe  suppression  is  about  30  dB. 

This  synthesis  program  outputs  the  normalized  current  weighting 
factors  (-l<a<+l)  and  does  not  determine  or  depend  on  the  technique 
being  used  to  achieve  the  weighting.  To  experimental ly  check  the 
synthesis  program,  four  sets  of  devices  were  constructed  (named 
CHECK5-A,  CHECK5-B,  CHECK7-A,  and  CHECK7-8)  which  are  shown  in 
Figure  3.7.  In  al 1  four  patterns,  the  transducers  were  3  mm  long 
and  30  um  wide  and  the  50  Ohm  transmission  lines  feeding  each 
element  were  equal  in  length  for  each  device.  This  was  done  to 
equalize  the  phase  for  each  bar.  The  devices  had  either  5  or  7 
bars  at  the  input  and  one  at  the  output.  MSW  delay  line  is  bilateral 
so  that  the  input  and  output  can  be  interchanged  when  the  direction 
of  the  magnetic  field  is  changed.  For  th»  two  cases  of  CHECK5-B 
and  CHECK7-B  the  patterns  were  designed  so  two  of  the  elements  could 
have  current  flowing  in  the  opposite  direction  of  the  current  in  the 
other  bars. 

In  addition,  different  attenuators  having  different  values  were 
connected  to  each  port  to  control  the  current  passing  through  each 
finger.  Using  equal  length  50  Ohm  coaxial  cables,  the  attenuators 
were  connected  to  an  8-way  p>ower  divider  where  its  unused  ports  were 
terminated  by  50  ohm  loads.  A  computer  program  (  Appendix  3)  was 
written  to  compute  the  insertion  loss  response  of  an  N  bar 
transducer  with  external  attenuators  for  current  weighting.  The 
theoretical  and  experimental  responses  for  these  devices  with 
different  cases  are  shown  in  Figure  3.8.  These  figures  demonstrate 
very  good  correlation  between  the  theory  and  experiment. 

The  next  step  was  to  construct  the  filter  and  use  a  current 
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CHECK7-A 


figure  3.7  See  text  for  explanations. 


Figure  3.  8  The  values  of  the  external  atternuators,  the  theoretical  predtiction 

of  the  frequency  response,  and  the  experimental  results  for 
several  devices  with  different  configurations.  (Notice  the  values 
of  the  attenuators  given  in  dB  for  each  case). 
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Figure  3.8  continued. 


Values  of  the  external  attenuators  (right),  theoretical  prediction  (top-left), 
and  the  experimental  results  (bottom-left). 


Figure  3.8  continued. 


In 


CEED 


CSZD 

■C 

<323 

<313 

■C2ZD 


Values  of  the  external  attenuators  (right), 
theortetical  predition  (top-left),  and  the 
external  results  (bottom-left). 

Figure  3.8  continued. 
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Values  of  the  external  attenuators  (top),  theoretical  prediction  (middle), 
and  the  experimental  results  (bottom). 

Figure  3.8  continued. 


Values  of  the  external  attenuators  (top),  theoretical  prediction  (middle) 
and  the  experimental  results  (bottom). 

Figure  3.8  continued. 
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weighting  technique  that  can  produce  any  current  weighting  factor  In 
the  range  of  - 1  to  + 1 . 

The  earl ler  method  of  weighting  the  transducer  by  thick  and 
thin  gold  plating  has  several  shortcomings.  First  of  al 1 ,  it  Is 
practically  limited  by  the  precision  of  controlling  the  thickness. 
Second,  one  can  only  achieve  positive  weighting  factors,  and  even 
then.  It  is  Impossible  to  have  a  practical  current  weighting  of 
zero. 

To  overcome  these  problems,  another  approach  to  the 
distribution  of  current  was  taken  by  using  an  open  gap  transducer, 
technique.  In  an  open  ended  transducer  bar  the  current  at  the 
position  of  the  gap  is  zero  and  Increases  In  the  sinusoidal  form  as 
distance  away  from  the  gap  increases.  If  the  length  of  the 
transducer  Is  very  small  In  comparison  with  the  wavelength,  the 
current  distribution  can  be  considered  to  be  linear.  If  a  transducer 
Is  fed  from  both  ends  and  a  gap  Is  placed  In  the  center  of  the 
transducer,  the  .net  current  of  the  transducer  will  be  zero.  By 
varying  the  position  of  the  gap  along  the  length  of  the  transducer. 
It  Is  possible  to  have  a  normalized  current  of  any  value  from  -1  to 
+1  (see  Figure  3.9).  Figure  3.10  shows  the  circuit  pattern  and  the 
frequency  response  of  a  21  bar  bandpass  filter  (named  FILTER-A) 
designed  to  have  100  MHz  bandwidth  centered  at  2.8  GHz.  Two  other 
designs  (FILTER-B  and  FILTER-C)  and  their  frequency  responses  are 
shown  In  Figures  3.11  and  3.12. 

All  these  filters  seemed  to  work  only  at  a  single  frequency  and 
when  the  tuning  was  attempted,  the  response  became  disturbed  as 
shown  In  Figure  3. 13-d.  The  following  explanation  Is  believed  to 
cause  the  distortion.  As  shown  In  Figure  3.13-a,  In  addition  to  the 
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SIN(X)/X  current  distribution  of  the  transducer,  there  is  also  an 
additional  current  distribution  due  to  the  transmission  line  feeding 
the  fingers.  The  combined  effects  are  shown  in  Figure  3.13-b,  and 
the  frequency  response  for  this  kind  of  current  distribution  is 
plotted  in  Figure  3.13-c.  This  is  very  similar  to  Figure  3. 13-d. 
Effects  of  the  addition  of  different  types  of  current  distributions 
on  the  main  SIN(X)/X  curve  are  shown  In  a  series  of  plots  in  Figure 
3.  14. 

A  power  divider  to  distribute  current  evenly  to  each  bar  was 
needed  to  overcome  this  problem.  A  simplified  current  divider  was 
used  to  design  a  16  finger  bandpass  filter  which  is  shown  in  Figure 
3.15.  The  tuning  capability  of  this  device  is  demonstrated  in 
Figure  3.16. 

Next,  to  suppress  the  other  harmonics  of  the  filter,  a  loop 
transducer  was  used  instead  of  the  single  bar  transducer.  As  shown 
in  Figure  3.17,  when  the  response  of  a  single  bar- loop  delay  line 
(curve  #2)  is  superimposed  on  the  frequency  response  of  a  weighted 
multibar-single  bar  filter  (curve  #  1),  the  result  is  suppression  of 
the  harmonics  (curve  #3).  The  procedure  to  theoretically  calculate 
the  1 oop-mu 1 t i bar  response  follows.  First,  program  1  of  Appendix  2 
is  run  for  single  bar-single  bar  mode.  Then  the  same  program  is  run 
for  loop- loop  mode.  The  two  outputs  are  then  used  by  program  4  of 
Appendix  2  to  obtain  the  linear  frequency  vs.  insertion  loss 
response  of  single  bar- loop  delay  line  with  zero  path.  The  last 
data  can  now  be  used  for  synthesizing  a  filter  in  the  manner 
described  earlier.  Figure  3.18  shows  the  circuit  pattern  and  a 
photograph  of  the  constructed  16  bar  filter  with  loop  transducer. 
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Figure  3.9  Current  Weighting  Filter 


68 


C.F  z  24  QHz 
A  A;  :  15  GHZ 


Figure  3.10  The  circuit  pattern  and  the  frequency 
response  of  FILTER-A. 
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Figure  3,1 1  The  circuit  pattern  (TOP),  theoretical  (MIDDLE)  and  the 
experimental  response  (BOTTOM)  of  the  device  FELTER-B.  The 
response  is  amplified  to  show  the  suppression  of  the  sidelobes. 
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Figure  3.12  The  circuit  pattern  (TOP),  theoretical  prediction  (middle)  and 
the  experimental  response  (bottom)  of  FDLTER-C.  The  output  is  ampliHed 
to  show  sidelobe  suppression. 
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Figure  3. 13  The  current  disctribudon  of  the  fingers  and  the  feeders 
(13-1)  combined  to  form  a  current  distribution  shown  in  figure  13-b 
The  predicted  frequency  response  (13-c)  is  very  similar  to  the 
experimental  response  of  the  devise  tuned  at  a  frequency  defferent 
from  the  one  designed  for. 
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Figure  3.14  The  theoretical  prediction  of  the  frequency  response  of 
a  00  MHz  bandpass  filter  when  the  feeders  contribute  t 
the  main  SIN(X)/X  current  distribution.  (3  different 
cooblnations  are  being  considered.) 


Figure  3. 14  continued. 
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pattern  of  a  16  bar  bandpass 
a  simplified  power  divider. 
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The  fre 
tuned  a 
500  Oe. 


X. 


Figure  3. 


17  Combining  Che  frequency  response  of  a  mulcibar-single  bar 
filter  with  a  single  bar-loop  delay  line  (TOP)  can  result 
in  suppression  of  the  harmonics. 
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Circuit  pattern  of  the  100  MHz  tunable  bandpass  filter. 


Photograph  of  the  constructed  16  element  bandpass  filter. 

Figure  3.18 
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The  frequency  response  and  its  tunability  are  shown  in  Figures  3.19 
and  3.20. 

To  improve  the  insertion  loss,  an  attempt  was  made  tc  narrow 
band  match  the  power  divider.  The  pattern  and  the  frequency 

response  is  shown  in  Figure  3.21.  As  can  be  seen  an  improvement  of 

about  10  dB  can  be  achieved  by  having  a  matched  power  divider. 
Another  filter  was  constructed  with  transducer  length  of  5  mm  and  was 
operated  at  higher  frequency  (4.5  -  5.5  GHz).  The  theoretical  and 
experimental  responses  are  shown  in  Figure  3.22.  With  this  device, 
the  minimum  insertion  loss  was  as  low  as  6  dB. 

A  matched  power  divider,  was  needed  to  gain  information  about 
the  transducer  coupling  at  different  frequencies  and  with  different 
lengths  of  transducer  being  under  the  YIG.  A  single  bar  I/O  delay 
line  was  used.  The  Sjj  response  of  the  device  at  2  GHz,  3  GHz,  4 
GHz  and  5  GHz  are  shown  in  Figure  3.23  as  points  A  through  D.  In 
this  case  the  YIG  was  on  the  entire  3  mm  length  of  the  transducer. 
Next,  the  effect  of  transducer  length  on  S j j  was  investigated  by 
using  25  urn  thick  YIG  covering  3  mm  (point  D),  2  mm  (point  E),  1  nnm 

(point  F)  and  0  mm  (point  G)  of  the  transducer  operating  at  5.0  GHz. 

As  shown  in  these  two  figures,  construct i ng  a  matched  power 
divider  is  a  difficult  task  because  of  the  complexity  of  the  load  on 
each  transducer. 

Finally,  instead  of  a  single  bar  in  output,  the  same  filter 
pattern  was  used  both  in  the  input  and  the  output  ports.  The 
theoretical  response  of  this  kind  of  filter  Is  compared  with  the  one 
which  has  single  bar  as  an  output  (Figure  3.24).  As  can  be  seen, 
theoretically  it  is  possible  to  suppress  the  side  lobe  as  much  as  45 
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10  ciB/dlv 


Theoretical  prediction  of  the  I.L.  vs  frequency. 


1.5  0.15  GHz/div  3.0 


The  experimental  I.L.  vs  frequency  response  of  the  device. 


Figure  3.19 
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2.0  0.1  GHz/div  3.5 

The  theoretical  simulation  of  the  tunability  of  the  filter. 
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The  experimental  results  of  the  filter  tuned  at  different  bias  fields. 

Figure  3.20 
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Figure  3.22  The  circuit  pattern  of  the  16  finger  filter  with  5  mm  transducer 

length  (TOP)  and  its  theoretical  simulation  of  frequency  response 
in  tunning  (BOTTOM-LEFT)  and  the  experimental  result 
(BOTTOM-RIGHT). 


Figure  3.23  TOP-  the  measured  S,,  of  a  single  bar  delay  line  with  the  YIG  covering 
entire  length  of  the  transducer  and  tuned  at  2  GHz  (point  A),  3  GHz 
(point  B),  4  GHz  (point  C)  and  5  GHz  (point  D). 

BOTTOM-  the  measured  S,,  response  of  the  same  device  tuned  at  5  GHz 
with  the  YIG  covering  entire  3  mm  transducer  (point  D),  covering 
2  mm  (point  E),  covering  1  mm  (point  F)  and  removed  (point  G). 
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dB.  But  when  a  device  was  made  (Figure  25)  with  a  semi matched  power 
divider,  the  results  were  not  satisfactory.  A  paper  describing  the 
theory  and  the  experimental  results  of  the  first  MSW  transversal 
filter  based  on  the  work  described  so  far,  was  presented  at  the 
lEEE-MTT  (1986)  meeting  and  was  published  In  the  proceedings  of  that 
meeting.  A  copy  of  this  paper  is  given  in  Appendix  4. 
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Figure  3.25  Pliocograph  of  the  device  with  synthesised  filter  pattern 
both  in  the  input  and  the  output,  and  the  frequency 
response . 
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^  .  DELAY  l_  I  NES 

4.1  S  Band  Linear  Delay  Line  The  research  for  the 
construction  of  linearly  dispersive  MSSW  delay  line  was  continued. 
Several  attempts  were  made  for  miniaturization  of  the  device  and 
improvement  of  its  insertion  loss  and  group  delay.  One  of  the 
devices  is  shown  in  Figure  4.1  along  with  the  measured  Sjj  and  S2J 
responses . 

The  magnet  package  shown  in  Figure  4.2a,  provided  by  Rockwell 
International  was  designed  to  house  the  MSSW  delay  line.  It  had  a 
constant  field  of  400  Oe  at  the  center.  The  device  (  Figure  4.2b) 
was  designed  to  have  the  transducers  located  in  the  center  of  the 
device.  Figure  4.2c  is  the  photograph  of  this  device  with  the 
graded  ground  plane  and  the  YIG  placed  on  it.  Figure  4.3  shows  the 
insertion  loss  vs.  frequency  response  in  the  frequency  range  of  2.4- 
3.9  GHz  (top- left),  the  reflection  coefficient  (Sjj)  for  2. 8-3. 2  GHz 
range  (top-right),  and  the  phase  -  magnitude  response  in  the  2.9- 
3.15  GHz  range  (bottom). 

Several  graded  ground  planes  were  tested  and  the  R.M.S. 
deviation  from  quadratic  phase  for  different  frequency  ranges  were 
calculated.  The  results  are  shown  in  Figure  4.4. 

Our  earlier  experience  had  shown  that  a  few  dB  improvement  in 
the  insertion  loss  can  be  achieved  if  the  graded  ground  plane  and 
the  ground  block  are  both  gold  plated,  and  if  indium  solder  is  used 
for  grounding  the  2^3  S'-'^strate  instead  of  silver  paste.  In  the 
final  stage,  when  the  repeated  measurements  of  the  insertion  loss 
response  and  the  group  delay  was  satisfactory,  the  graded  ground 
plane  was  permanently  soldered  to  the  block  and  the  YIG  strip  was 
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•••••MSSW5***** 

B.W=  18000  TO  3.0550  GHz 
A»  181.4271  8=^70.4292 
Phue(QUADRATIO.JA(F2-F^)-B*(F-F0) 

THE  R.M.S.  DEVIATION  FROM  QUADRATIC  PHASE  IS;  20.8681 


•••♦•MSSW5***** 

B.W=  18500  TO  3.1850  GHz 
A=  1541979  B=-390.4362 
PhMe(QUADRATI0.5A(F2-F02)-B(F-Fo) 

THE  R.M.S.  DEVIATION  FROM  QUADRATIC  PHASE  IS:  7.5878 


B.W=  19000  TO  3.1550  GHz 
A=  133.6094  B=-328.6460 
Ph«se<QUADRATIC)=.5A(F2.IH)2)-B(F-Fo) 

THE  R.M.S.  DEVIATION  FROM  QUADRATIC  PHASE  IS:  13.8825 


B.W-  19450  TO  3.2000  GHz 
A-  1123049  B— 264.0575 
Phue(QUADRATI03A(F2.F^).B(F-Fo) 

THE  R3I.S.  DEVIATION  FROM  QUADRATIC  PHASE  IS:  16.0892 


Rgnre44  continued. 
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secured  with  glue  applied  at  both  ends.  This  linearly  dispersive 
(positive  slope)  MSSW  delay  line  was  Intended  to  be  cascaded  with 
another  linearly  dispersive  (negative  slope)  MSBVW  delay  line  to 
produce  an  electronically  tunable  constant  delay  line.  The  MSBVW 
device  was  constructed  at  Rockwell  International.  A  paper 
describing  the  theory  of  operation  and  the  experimental  results  was 
presented  at  the  Phased  Arrays,  1985  Symposium  and  was  also 
published  In  the  1986  March  Issue  of  the  MICROWAVE  JOURNAL.  A 
reproduction  of  this  paper  Is  given  In  Appendix  5. 

4.2  Ka  Band  Delay  Lines  To  explore  the  usefulness  of  MSW 
devices  at  higher  frequencies,  the  response  of  simple  delay  line  at 
Ka  band  (26.5-40  GHz)  was  Investigated.  The  first  sample  was  an 
open  ended  single  bar  I/O  delay  line  with  transducer  length  of  3  mm 
and  10  urn  width  constructed  on  quartz  substrate. 

Figure  4.5  shows  photographs  of  the  S2J  for  the  zero  line  (no 
device)  and  also  the  direct  break-through  (device  with  no  H  field). 
Figure  4.6  Is  the  frequency  response  of  the  device  operated  In  the 
forward  volume  wave  mode  for  different  YIG  widths  and  different 
positions.  The  same  procedures  are  repeated  for  the  backward  volume 
wave  mode  and  the  results  are  shown  In  Figure  4.7.  All  the  pictures 
In  Figures  4.5,  4.6  and  4.7  were  taken  for  center  frequency  of  27 
GHz  and  frequency  sweep  of  1  GHz  with  50  um  thick  YIG. 

The  tunablllty  of  the  device  In  the  frequency  range  of  27  GHz- 
33  GHz  Is  shown  Is  the  series  of  photographs  displaying  Sjj  and  Sj2 
In  Figure  4.8.  For  these  pictures  the  1.5  mm  wide  and  50  um  thick 
YIG  covered  the  second  half  (close  to  the  transmission  line)  of  the 
3  mm  long  tra.srlticer .  For  a!!  the  p’ct'j'“es  the  frequency  sweep  was  I 
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Figure  4. 


CF  27  GHz 
A  F  1  " 


S'l  response  fur  ,':or('  li.io  (no  dovico)  and  direcc- 
( device  without  the  field). 
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CF  2  7  OHZ 


Figure  4.6  I.L  vs  frequency  response  for  different  YIG  widths  and  positions 
for  MSFVW  mode. 
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GHz.  The  second  device  was  a  150  um  loop  I/O  delay  line  (with  1 
mm  center/center  spacing)  constructed  on  quartz  substrate.  A 
photograph  showing  the  pattern  and  also  the  frequency  response  of 
the  device  Is  shown  In  Figure  4.9.  In  the  picture  showing  the 
response,  the  top  line  Is  the  zero  line,  the  bottom  line  Is  the 
direct  break-through  and  the  middle  line  Is  S2 j  of  the  device,  all 
in  the  frequency  range  of  27.5-28.5  GHz  operated  In  the  MSFVW  mode 
with  50  urn  thick  YIG. 

After  subtracting  the  zero  line  from  the  j  of  the  MSFVW  mode 
In  the  frequency  range  of  27-29  GHz,  a  sample  plot  of  Insertion  loss 
vs.  frequency  for  another  device  is  shown  In  Figure  4.10.  This 
figure  shows  there  is  a  promise  of  designing  MSW  delay  lines  at  high 
frequenc 1 es . 

Another  device  (Figure  4.11)  was  designed  with  open  ended 
transducers  of  width  30  um  and  length  4  mm.  In  order  to  reduce  the 
mismatch,  special  semi-rigid  cables  with  type  K  connectors  were  used 
to  Interface  the  device  with  the  test  instrument. 

Several  experiments  were  performed  with  YIG  strips  of  different 
thicknesses  and  widths,  and  positioned  at  various  places  on  the 
transducers.  The  results  showed  that  narrower  strips  (less  than 
I  mm)  with  YIG  thicknesses  >  50um  positioned  2/3  of  the  way  away 
from  the  open  end  of  the  transducer  resulted  in  the  lowest  insertion 
loss  and  the  smoothest  response.  Figure  4.12  shows  the  insertion 
loss  vs.  frequency  and  phase  response  operated  in  the  MSFVW  mode. 
The  pictures  are  for  the  center  frequency  of  30  GHz  and  a  bandwidth 
of  1.5  GHz.  For  the  same  center  frequency,  but  with  the  sweep 
bandwidth  of  500  MHz,  the  magnitude  and  the  phase  responses  in  the 
same  mode  are  shown  in  Figure  4.13.  The  magnitude  and  the  phase 


101 


responses 
of  MSSW 


of  the  dev 
are  shown  I n 


ce  operated  in  the  forward  and  reverse 
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UMMARY 


In  summary  the  primary  accomp I f shments  of  this  final  program 
period  have  been: 

1.  A  variable  time  delay  device  with  an  R.M.S.  phase  error  of 
ab  low  as  6.40°R.M.S.  has  been  built  and  tested.  Time  variation  of 
42  nsec  was  achieved.  the  design  procecures  for  dispersion  control 
utilizing  a  variable  and  backward  wave  delay  line  with  linearized 
dispersion,  and  thus  matched  dispersive  line  have  been  developed 
with  significant  potential  in  signal  processing  applications. 

2.  The  current  weighting  technique  using  open  gap  transducers 
to  realize  MSSW  filters,  has  demonstrated  good  correlation  between 
the  theoretical  prediction  and  measured  response.  Using  this 
technique,  filters  were  realized.  Tunable  bandpass  filters  with 
side  lobe  suppression  greater  than  20  db  and  passband  ripple  of  less 
than  3  db,  and  insertion  loss  of  6  db  have  been  realized. 

3.  Ka  Band  delay  MSFVW  lines  with  insertion  losses  of  less 
than  16  db  and  bandwidth  of  greater  than  1  GHz  have  been  realized. 
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Figure  4.10  Plot  of  I.L  vs  F  for  device  operating  in  the  MSFVW  mode 
corrected  for  the  zero  line. 
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Figurr'  4.13  The  expanded  view  of  Figure  4.1.2. 
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Insertion  Loss,  IL^  (dB) 


Frequency  £  (GHz) 


The  frequency  response  for  a  zero  path  single  bar  I/O 
transducer,  A(ui),  and  the  desired  response  of  a  multlbar 
bandpass  fitter,  H(w). 


G(a))  -  "  ^n-1  )e  (1) 


Gifii)  aarray  response 
N  •number  of  transducer  elements. 

ip^Cu))  "nth  element  response. 

a^  •weighting  factor  for  the  nth  element. 

kCim)  •transmission  response  of  a  single-bar /single-bar 
.  zero-path  delay  line. 

K(u)  •  a+j  B  (a  and  Bare  loss  factor  and  wave  number, 
respectively) . 

Z  •path  length  for  the  nth  element. 

H(u))  •desired  frequency  response. 
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APPENDIX  2 

"LIST  OF  THE  PROGRAMS  FOR  SYNTHESIZING  BANDPASS  FILTER" 


116 


o  o  o 


PROGRAM  //I 


C  MSSW  DISPERSION  RELATION  CALCULATES  152  U-K  POINTS  AND  305  W-K  ' 

C  POINTS  FOR  HILBERT  TRANSFORM  INTEGRATION 

C  MSSH  DISPERSION  RELATION,  GROUP  VELOCITY  AND  TIME  DELAY, PHASE  VEL. 
C  DESIGNED  FOR  STRUCTURE  A/  B,  REVERSE  (+/-)  DIRECTION 

^«***«****««*«««*«**********«««k««*«««*«r«»*««**«******«*)i* 

C  EDITTED  AND  MODEFIED  ON  7/24/84.  Y.J.  ATAIIYAN 

DIMENSION  FRQ(152) ,RK1(152) ,RK2(152) ,TM1(152> ,TM2(152) ,QF(305) ,QK1 
1(305), QK2(305) 

DIMENSION  RI1(152),RI2(152),RIC(152),XI1(152),XI2(152), 
1XIL(152),T01(20) 

DIMENSION  R1(152),R2(152),RMX1(152),RMX2(152),XMX1(152), 
1XMX2(152),QRM(305),QRN(305),FQN<152),WAVL(152),ILMA(152),ILMB(152) 
DIMENSION  FQ(152),TD1(152),TD2{152),XILMA(152) 

REAL  ILMA.ILMB 

REALMS  VP1,VP2,RKF,RKR,TM1,TM2 

REALMS  FRQ , R1 , Rl , QF , QRl , QK2 , FRQ3 , FRQ4 , FRQl , FRQ2 , DF , RFRQ , QRM , F , 2F , R 
IKl , RK2 , QRN , RMXl , RMX2 , ARGUE , HI , H2 ,H3 , SUMl , SUM2 , SUM , SIMP2 , RM , RMl , RM2 
1 , RND . OH , RM , RML2 , XML2 , GM , XMXl , XMX2 
REAL*8  FQ , R22 , TD2 , XMl , XM2 , XM , EK , 
1ZC,BETA.AC,AR.AR1,AR2,AX,AX1,AX2.CZC,CBT,CAC,CAR, 
2CXM,ALPH,AL,BL,RIN,XIM,RI1,RI2,RIC,XI1,XI2,XIL,PT1,PT2,CL1,CL2 
3,XL1,XL2,X,CST 

REAL*8  T,0W,RNMAX,PN,U1,U2,AU1,AU2, 

1RK,RKS,RD,RT,BT,FT,DELT,EXT,A,B,C,D,PC,P1,P3,P21,P22,P2,P, 
,DET1,DET2,E,G,RTT,RTU,TT,R11 
REAL*8  SSECH 

REAL*8  DELTF,A1,A2,RTL,RT1,TL, 

,T1,DL,D1,FL,F1,TX,DX,FX,V0,V1,V7,V3,VG1,VG2,TDD,TDT 
COMMON  F,RKF,RKR,RM1,RM2 
COMMON  DO,TO,BO,RLGH.TU,SO,UO 
COMMON  HI,GA,JCONF,NMAX,FPM 
DSECH (X ) »DSqRT( 1 . -DTANH ( X ) *DTANH ( X ) ) 

(;  ***«»  DEFINITION  OF  THE  TERMS  ************************** 

C  DO=YIG  THICKNESS  (METERS) 

C  TO*  AL203  THICKNESS  (") 

C  BO*  TRANSDUCER  HIDTH  (") 

C  RLGH*  TRANSDUCER  LENGTH  (") 

C  TU*  YIG/TRANSDUCER  AIR  GAP  (") 

C  SO*  SEPARATION  BETWEEN  TRANSDUCERS  (") 

C  UO*  ?!!I 
C  HI*  BIAS  FIELD  (OE) 

C  GA*  GYR0HA6NETIC  RATIO  (HZ/OE) 

C  Z0«  CHARACTRISTIC  IMP.  OF  THE  TRANSDUCER. 

SO*WAVELENGTH  AT  CENTER  FREQ.  IT  DOES  NOT  HAVE  ANY  EFFECT 
ON  SINGLE  BAR,  FOR  LOOP  TRANSD.  SPACINGS  ARE  SO/2,  FOR 
PARALLELL  IT  IS  EQUAL  TO  SO. 

C  CONDITION  JUMP  ON  TRANSDUCER  CONFIGURATIONS  (  1  FOR  SINGLE  BAR,  2 
C  FOR  PARALLEL  ODD  BARS,  3  FOR  PARALLEL  EVEN  BARS,  4  FOR  MEANDER 

C  STRIPS,  5  FOR  PI  0.5  CENTER  SPACING,  6  FOR  PI  1.5  CENTER  SPACING) 

C  NMAX  IS  THE  STRIP  NUMBER  FOR  MSU  TRANSDUCERS 
C  INPUT  THE  FILE  NAME  FOR  STORING  THE  DATA 
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OPEN (UNIT»23, DEVI CEa'DSK', DIALOG) 

TYPE  1200 

1200  FORHAT( IX, 'ENTER  YIG  AND  ALUMINA  THICK  IN  MICRON') 
ACCEPT  1201, DO, TO 

1201  FORMAT  (2F) 

D0»D0*1.0E-6 

T0»T0*1.0E-6 

FPM=1750. 

TYPE  1202 

1202  FORMAT (IX, 'ENTER  EXT  H  FIELD') 

ACCEPT  1203, HI 

1203  FORMAT(IF) 

GAs2.8E6 

GM»GA»FPM 

TYPE  1204 

1204  F0RHAT(1X,'TYPE  UIDTH  AND  LENGTH  OF  TRANSD.  IN  MICRON') 
ACCEPT  1201,B0,RLGH 

TYPE  1210 

1210  FORMAT (IX, 'ENTER  I/O  TRANSD.  SPACINGS  IN  CM') 

ACCEPT  1203, PATH 
PATH-PATH* l.OE-2 
TYPE  1205 

1205  FORHAT( IX, 'ENTER  AIR  GAP(HICRON),  AND  LAMBDA  AT  CENTER' 
•,'FREQ.  (MICRON)') 

ACCEPT  1201,TU,S0 
BO-BO* l.OE-6 
RLGH«RLGH*1.0E>6 
TU-TU*1.0E-6 
SO-SO* l.OE-6 
U0-1.256637E-6 
TYPE  1206 

1206  FORMATdX,'  ENTER  CHARAC.  IMP.  OF  TRANSD  (ZC)') 

ACCEPT  1203,ZC 

TYPE  1207 

1207  FORMATdX, 'ENTER  COND  CODE  AND  *  OF  TRANSDS  (I  I)') 

TYPE  1200 

1208  FORMATdX,'  CC  FOR  SINGLE-1,  PARAL  OOD-2,  PARAL  EVEN-3 
*,MEAND-4,  .5  SPACING  PI-5,  1.5  SPACING  PI-6') 

ACCEPT  1209,JC0NF,NMAX 

1209  F0RHAT(2I) 

TYPE  1213 

1213  FORMATdX, 'IL  IN  OB  (TYPE  1)  OR  LINEAR  IL  (TYPE  0)?') 
ACCEPT  1214,NYJA 

1214  FORHAT(I) 

OH-HI/FPM 

FRQ3-DSQRT(0H*(0H-fl. )  )*GM 
FRQ4-(0H-M).5)*GN 
FRQ1«FRQ3*1.01 
FRQ2»FRQ4*0.99 

C  FREQUENCY  SAMPLING  305  POINTS  FOR  HILBERT  TRANSFORM 
qF(l)«FRq3 
qF(3)-FRqi 

qF(2)-0.5*(FRq3-»-FRqi) 

qF(103)-0.5*(FRq2+FRqi) 

qF(305)-FRq4 
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QF(303)sFRQ2 

QF ( 304 )«0 . 5» ( FRQ2+FRQ4 ) 

DF=(QF(103)-QF(3))/100.0 

DO  100  I«l,99 

RI»I 

J«I+3 

100  QF(J)»QF(3)+RI*DF 
DO  101  I»l,199 
RI=I 

J>I-t-103 

101  QF<J)»QF(103)+RI»DF»0.5 

C  FREQUENCY  SAMPLING  152  POINTS  FOR  DISPERSION  RELATION 
FRQ ( 1) =FRQ  3+0 . 7  5* ( FRQ 1 -FRQ3 ) 

FRQ(2)=FRQ1+0.5*DF 
DO  102  Isl.49 
RI=I 
J»I+2 

102  FRQ(J)3FRQ<2)+RI*DF*2.0 
FRQ(52)=QF(103)+DF*0.25 
DO  103  1=1,99 

RI=I 

J=52+I 

103  FRQ(J)=FRQ(52)+RI*DF 
FRQ(152)=FRQ2+0.25»(FRQ4-FRQ2) 

C  DISPERSION  RELATION 

DO  130  H*1,1S2 
F*FRQ(M) 

OW-F/GM 

Ul=l .0-OH/(OW*OW-OH»OH) 

U2=0W/(0W*0W-0H*0H) 

C  S=+l 

AUU1+U2+1. 

A2=U1-U2-1 . 

C  REGI  FALSI  ITERATION  TO  FIND  INITIAL  GUESS  FOR  NEVTON  ITERATION 
RTL=0.1*T0 
RT1=1.0E5»T0 
DO  20  1=1,500 

IF((I-500).EQ.0)  WRITE(23,13) 

TL=DTANH(RTL) 

T1=0TANH(RT1) 

DL=DEXP ( -2 . *DO*RTL/TO ) 

D 1-DEXP ( -2 . *D0*RT1 /TO ) 

FL»(U1-U2+TL) *A1-(UI+U2-TL ) »A2»DL 
F1«(UI-U2+T1 )»A1-(U1+U2-T1 )»A2*D1 
X=(RTL*F1-RT1*FL)/(F1-FL) 

TX=OTANH(Z) 

DX=DEXP(-2.*D0*X/T0) 

FX»(U1-U2+TX)*A1-(U1+U2-TX)*A2*DX 
IF(DABS<FX)-1.0B-7)  51,51,31 
31  CLX=FX»FL 
CUX=FX*F1 

IF(CUX.CT.O.)  RT1»X 
IF(CLX.GT.O.)  RTL»X 
20  CONTINUE 

C  NEVTON  ITERATION 
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51  DO  90  L=l,500 

IF((L-500).EQ.0)  WRITE(23,19) 

TX«OTANH(Z) 

DX«DEXP(-2.*D0*X/T0) 

FX»(U1-U2+TX)*A1-(UI+U2-TX)»A2*DX 
DF=DX»A2*(+2.*D0*(U1+U2-TX)/T0+0SECH<X)*»2)+A1*DSECH(X)**2 
IF(DABS(FX)-1.0E-10)  92,92,91 

91  X«X-FX/DF 
90  CONTINUE 

92  R11=X/T0 

C  PHASE  VELOCITY,  GROUP  VELOCITY  AND  DELAY  TIME  CALCULATION 
RK1(H)-R11 
FN1»FX 
TX>DTANH(X) 

DX=DEXP(+2.»D0*X/T0) 

V0«3. 078760801 D10*T0 

V1=DX*AI*(+2.*DO*(U1-U2+TX)/TO+DSECH(X)»*2)+A2*DSECH(X>*»2 

V2»(<U1+U2-TX)-A1*DX)/((0H-0H)**2) 

V3«(A2-(U1-U2+TX)*DX)/((0H+0H)**2) 

VG1»V0*V1/<V2-V3) 

TDD=PATH/VG1 

TM1(B)=TDD 

VPl=6.2a3185307*F/Rll 
C  Sa-1 

A1»U1+U2-1 . 

A2*U1-U2+1 . 

C  REGI  FALSI  ITERATION  TO  FIND  INITIAL  GUESS  FOR  NEWTON  ITERATION 
RTL«0.1»T0 
RT1=1.0E5»T0 
DO  209  I»l,500 

IF((I-500).EQ.0)  WRITE(23,139) 

TL=DTANH(RTL) 

T1=DTAHH(RT1) 

DL=DEXP ( -2 . »DO*RTL/TO ) 

D  UDEXP  <  -2 .  *D0*RT1  /TO ) 

FLxDL* (U1-U2-TL ) »A1-(U1+U2+TL ) *A2 
F1=D1*<U1-U2-T1)»A1-<U1+U2+T1)*A2 
X»(RTL*F1-RT1*FL)/(F1-FL) 

TX>DTANH(X} 

DX-DEXP ( -2 . »D0*X/T0 ) 

FX»DX* (U1-U2-TX) *Al-<  U1+U2+TX) * A2 
IF(0ABS(FX)-1.0E-7)  519,519,319 
319  CLX»FX*FL 
CUX«FX»F1 

IF(CUX.GT.O.)  RTl-X 
IF(CLX.GT.O.)  RTL-X 
209  CONTINUE 
C  NEWTON  ITERATION 
519  DO  909  L»l,500 

IF((L-500).EQ.0)  WRITE(23,199) 

TX«DTAMH(X) 

DX«DEXP(-2.*D0*X/T0) 

FX»DX*(U1-U2-TX)*A1-(U1+U2+TX)*A2 

DF«DX»A1*(-2.»00*(U1-U2-TX)/T0-DSECH<X)»*2)-A2*DSECH<X)»»2 
IF(0ABS(FX)-1.0E-10)  929,929,919 
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919  X-X-FX/DF 
909  CONTINUE 
929  R22=X/TO 

C  PHASE  VELOCITY,  GROUP  VELOCITY  AND  DELAY  TIME  CALCULATION 
RK2(H)aR22 
FN2sFX 
TX=DTANH(X) 

DX=DEXP(-2.*D0*X/T0) 

VO>3.078760801D10*TO 

VI«DX*A1*(-2.*D0*(U1-U2-TX)/T0-DSECH(X)**2)-A2*DSECH{X)**2 
V2»( (Ul+U2+TX)-Al*DX)/( (0H-0H)»*2) 

V3» ( A2-( U1 -U2-TX ) *DX ) / ( ( OW+OH ) **2 ) 

VG2=V0*V1/(V2-V3) 

TDT=1.0E-2/VG2 

TM2(M)»TDT 

VP2=6.283185307*F/R22 
130  CONTINUE 

C  HILBERT  TRANSFORM  H-K  POINTS 

DO  5130  M-2.304 
FxQF(M) 

OH»F/GM 

Ul»l .0-OH/(OW»OW-OH*OH) 

U2=0W/(0U*0W-0H*0H) 

C  S=+l 

AlsUl-t>U2-fl. 

A2=U1-U2-1 . 

C  REGI  FALSI  ITERATION  TO  FIND  INITIAL  GUESS  FOR  NEWTON  ITERATION 

RTLa0.l*T0 
RTl»1.0E5n0 
DO  520  I»l,500 

IF((I-500).EQ.0)  WRITE(23,13) 

TL=DTANH(RTL) 

T1=DTANH(RT1) 

DL=DEXP(-2.*D0*RTL/T0) 

DUDEXP(-2.»D0*RT1/T0) 

FL= ( U1-U2+TL ) »A1-(U1 +U2-TL ) * A2*DL 
F1«(U1-U2+T1)»A1-<U1+U2-T1)*A2*01 
X=(RTL*F1-RT1»FL)/(F1-FL) 

TXsDTANHCX) 

DX=DEXP(-2.*D0»X/T0) 

FX»(U1-U2+TX)*A1-(U1+U2-TX)*A2*DX 
IF(DABS(FX)-1.0E-7)  551,551,531 
531  CLXaFX»FL 
CUX»FX*F1 

IF(CUX.GT.O.)  RT1»X 
IF(CLX.GT.O.>  RTL«X 
520  CONTINUE 
C  NEWTON  ITERATION 
551  DO  590  L>1,500 

IF((L-500).EQ.0)  WRITE(23,19) 

TX«0TANH(X) 

DX«DEXP(-2.»D0*X/T0) 

FX«(U1-U2+TX)»A1-(U1+U2-TX)*A2*DX 
DF»DX*A2*(+2.*D0*(U1+U2-TX)/T0+DSECH(X)»*2)+A1*DSECH(X)»*2 
IF(DABS(FX)-1.0E>10)  592,592,591 
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591  X»X-FX/DF 
590  COMTIMUE 

592  R11=X/T0 
QKl(n)»Rll 
FMUFX 

C  S»-l 

AlaUl+U2-l. 

A2=U1-U2+1 . 

C  REGI  FALSI  ITERATION  TO  FIND  INITIAL  GUESS  FOR  NEWTON  ITERATION 
RTL=0.1*T0 
RT1»1.0E5*T0 
DO  5209  1=1,500 
IF((I-500).EQ.0)  URITE(23,139) 

TL=DTANH(RTL) 

T1=DTANH(RT1) 

DL=DEXP(-2.*D0*RTL/T0) 

D1=DEXP(-2.*D0»RT1/T0) 

FL=DL*(U1-U2-TL)*A1-(U1+U2+TL)*A2 

F1=D1*(U1-U2-T1)*A1-(U1+U2+T1)*A2 

X=(RTL*F1-RT1»FL)/(F1-FL) 

TX=DTANH<X) 

DX=DEXP(-2.*D0*X/T0) 

FX=DX*(U1-U2-TX)»A1-(U1+U2+TX)*A2 
IF(DABS(FX)-1.0E-7)  5519,5519,5319 
5319  -CLX=FX*FL 
CUX=FX»F1 

IF(CUX.GT.O.)  RT1»X 
IF(CLX.GT.O.)  RTL=X 
5209  CONTINUE 
C  NEWTON  ITERATION 
5519  DO  5909  L=l,500 

IF((L-500).EQ.O)  WRITE(23,199) 

TX=DTANH<X) 

DX=DEXP<-2.*D0*X/T0) 

FX=DX»(U1-U2-TX)*A1-(U1+U2+TX)*A2 

DF=DX*A1*(-2.*D0*(U1-U2-TX)/T0-DSECH(X)**2)-A2*DSECH(X)*»2 
IF(DABS(FX)-1.0E-10)  5929,5929,5919 
5919  X=X-FX/DF 
5909  CONTINUE 
5929  R22=X/T0 
qR2(H)=R22 
FN2»FX 

5130  CONTINUE 

C  nssw  RADIATION  IMPEDANCE  RH  AND  XH/  RH  AS  A  SUBPROGRAM,  XM 

C  CALCULATED  FROM  HILBERT  TRANSFORM  OF  XM  BY  USING  2ND  ORDER  CLOSED 

C  FORM  SIMPSON  RULE  FOR  INTEGRATION 
GM-GA»FPM 
OH»HI/FPM 

FRQ3=DSQRT(0H*(0H+1 . ) )*GM 
FRQ4=(OH+0.5)*GM 
FRqi=FRQ3*1.01 
FRq2»FRq4*0.99 
DF=<FRq2-FRqi)/200.0 
RFRq»FRq4/FRq3 
qF(l)»i.o 
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QF(305)»RFRQ 

QRH(1)»0. 

QRM(305)»0. 

QRM(1)=0. 

QRN(305)30. 

DO  200  1=2,304 
F»QF(I) 

RKF=QR1<I) 

RKR=QK2(I) 

CALL  P6RH 
QRH(I)=RH1 
QRN(I)=RH2 
QF(I)=QF(I)/FRQ3 
200  CONTINUE 

DO  400  1=1,152 
F=FRQ(I) 

ZF=F/FRQ3 

RKF=RK1(I) 

RKR=RK2(I) 

CALL  PGRH 

RnXl(I)=RHl 

RHX2(I)=RH2 

ARGUE= ( ZF+1 . ) * ( RFRQ-ZF ) / ( ( ZF-1 . ) * ( RFRQ+ZF ) ) 

RH«Rni 

RN=RH2 

C  2ND  ORDER  CLOSED  FORM  OF  SIMPSON  RULE  FOR  INTEGRATION 
Hl= < QRM ( 1 ) -RM) / ( QF ( 1) *QF ( 1 ) -2F*ZF ) 
H2=<QRM<2)-RM)/(QF(2)*QF(2)-ZF*2F) 

H3=( QRM<  3 )-RM) / ( OF ( 3 ) *QF ( 3 ) -ZF*ZF ) 

SUM1=( FRQ1-FRQ3 ) • (Hl+4 . *H2+H3 ) / (6 . *FRQ3  > 
Hl=(qRM(303)-RM)/<QF(303)*QF(303)-ZF*ZF) 
H2=(QRn(304)-RM)/(QF(3C4)*QF(304)-ZF*ZF) 
H3=(QRM<305)-RM)/(QF(305)*QF(305)-ZF*ZF) 

SUM2=( FRQ4-FRQ2 ) * (Hl+4 . *H2+H3 ) / ( 6 . *FRQ3 ) 

SUM=SUM1+SUM2 
DO  401  K=l,50 
JK=2*<K+1) 

J1=JK-1 

J2=JK+1 

H1=(QRM(J1)-RM)/(QF(J1)»QF(J1)-ZF*ZF) 

H2= ( QRM ( JK ) -RM ) / ( QF ( JK ) *qF  <  JK ) -ZF*ZF ) 
H3»(qRM(J2)-RM)/(qF(J2)*qF(J2)-ZF*ZF) 

SIMP2=DF*(Hl+4 . »H2+H3 ) / ( 3 . *FRq3 ) 

401  SUM=SUM+SIMP2 
DO  402  K«51,150 
JK»2*(K+1) 

J1»JK-1 

J2=JK+1 

Hl»(qRM(Jl)-RM)/(qF(Jl)*qF<Jl)-ZF»ZF) 

H2= ( qRM( JK ) -RM ) / ( OF ( JK) *qF ( JK) -ZF»ZF ) 
H3=(qRH(J2)-RM)/(qF(J2)*qF<J2)-ZF*ZF) 

SIMP2*DF* (Hl+4 . *H2+H3 ) *0 . 5/ ( 3 . »FRq3 ) 

402  SUM«SUn-l'SIMP2 

XMXl(I)=0.318098862»(RM*DL0G(ARGUE)-*-2.*ZF*SUM) 

C  2ND  ORDER  CLOSED  FORM  OF  SIMPSON  RULE  FOR  INTEGRATION 
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Hla(QRN(l)-RN)/(QF(l)*QF(l)-ZF»ZF) 

H2=(QRN(2)-RN)/(QF(2)*QF(2)-ZF»ZF) 

H3»(QRN(3)-Rll)/(qF(3)*QF(3)-ZF*ZF) 

SUM1*( FRQ1-FRQ3 ) * (Hl+4 . »H2+H3 ) / ( 6 . *FRQ3 ) 

Hl=(  QRJI(  303  )-RN )  /  ( QF  ( 303 )  *QF  ( 303  )-ZF*ZF  > 
H2a(QRN(304)-RH)/(QF(304)*QF(304)-ZF*ZF) 
H3a(QRN(305)-RN)/(QF(305)*QF(305)-ZF*ZF) 

SUM2=( FRQ4-FRQ2 ) * (Hl+4 . *H2+H3 ) / ( 6 . *FRQ3 ) 

SUMxSUMl+SUM2 
DO  550  K«1.50 
JK=2*(K+1) 

JUJK-l 

J2*JK+1 

Hlx(QRJI(Jl)-RM)/(QF(Jl)*QF(Jl)-ZF*ZF) 

H2a(QRN(JK)-RM)/(QF(JK)*QF(JK)-2F*ZF) 

H3=(QRN(J2)-RH)/(QF(J2)*QF(J2)-ZF*ZF) 

SIMP2=DF* (Hl+4 . *H2+H3 ) / ( 3 . *FRq3 ) 

550  SUK=SUM+SIMP2 
DO  555  Rx51,150 
JKa2*(K+l) 

J1=JK-1 

J2=JK+1 

Hlx(qRH<Jl)-RM)/(qF(Jl)*qF(Jl)-ZF*ZF) 

H2=(qRN(JK)-RN)/<qF(JK)*qF(JK)-ZF*ZF) 

H3»(qRH<J2)-RM)/(qF(J2)*qF(J2)-ZF*2F) 

SIMP2=DF* (Hl+4 . *H2+H3 )*0 .5/(3. *FRq3 ) 

555  SUMaSUM+SIMP2 

Xia2  ( I)xO .  318098862*  ( RN*DL0G  (ARGUE  )+2 .  *ZF*SUll) 

400  COMTINUE 

C  nSSV  FREqUEMCY  RESPONSE  PROGRAM  TO  CALCULATE  TRANSMISSION 

C  LINE  PARAMETERS,  CIRCUIT  ELEMENTS.  INPUT  IMPEDANCE,  AND  INSERTION 

C  LOSS 

RNxNMAX 

ND=NMAX/2 

RNDxND 

CST=1. 

DO  105  J=l,152 
F=FRq(J) 

TM9«TM1(J) 

TM8=TM2(J) 

RMlxRHXKJ) 

RM2xRMX2(J) 

XMlxXMXKJ) 

XM2-XMX2(J) 

RM>RM1-«'RM2 

XMxXHl>XH2 

BETA>3.636102608D-8«F 

ACa6.32D-7*0SqRT(F)/(B0*ZC) 

RML2«RM*RLGH*0.5 

XML2»XM*RLGH*0.5 

AR»0.5*RN/ZC 

ARlx0.5*RMl/ZC 

AR2«0.5*RM2/ZC 

AX«0.5*XM/ZC 

AXlx0.5*XMl/ZC 
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AX2=0.5*XM2/ZC 

GO  TO  (81,82,82,93,94,94),JCOMF 

81  CZC=l. 

CBT*1 . 

CACsl . 

CARsl. 

CXM=1 . 

GO  TO  99 

82  CZC=l./RN 
CBT=1 . 

CAC=1 . 

CAR»1./RM 
CXM=1./RN 
GO  TO  99 

93  CZC=1. 

CBT=RM 
CAC=RN 
CAR»1./RN 
CXM»1./RN 
GO  TO  99 

94  CZC=l./RND 
CBT»2. 

CAC=2. 

CARsl./RN 

CXM=1./RN 

C  TRANSMISSION  LINE  PARAMETER  MODIFICATIONS  FOR  TRANSDUCER  ARRAYS 
99  ZC=C2C*ZC 
BETAaCBT*BETA 
ACaCAC*AC 
AR=CAR*AR 
AR1=CAR»AR1 
AR2»CAR*AR2 
AX=CXM*AX 
AX1=CXM*AX1 
AX2=CXM*AX2 
ALPHsAR-hAC 
BET-BETA+AX 
AL=2.»ALPH*RLGH 
BL=2.0»BET»RLGH 

C  INPUT  IMPEDANCE  AND  CIRCUIT  ELEMENTS 

RIN>ZC*DTANH ( AL)  / ( 1 . -t-DCOS ( BL ) *DSECH( AL } ) 

XIN=ZC*DSIN(BL)*DSECH(AL)/(l.+DCOS(BL)»DSECH(AL)) 

RIl(J)aRIM*ABl/ALPH 

RI2(J)aRIN*AR2/ALPH 

RIC(J)aRIN*AC/ALPH 

XIl(J)aXIN*AXl/BET 

XI2(J)aXIN*AX2/BET 

XIL(J)»XIN»BETA/BET 

C  POWER  RATIO,  CONVERSION  LOSS  AND  INSERTION  LOSS 
PTla200.0*RIl(J)/((RIN+50.0)»*2.+XIN**2) 
PT2a200.0*RI2(J)/({RIN+50.0)»*2.+XIN**2) 

CLUIO.*DLOG10(PT1) 

CL2alO.*DLOGlO(PT2) 

XL1»2 . *CLl-76 . 4E+6*TM9*CST* 1 . 5 
XL2a2. •CL2-76. 4E+6*TM8*CST*1 . 5 
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IF(MYJA.EQ.l)  GO  TO  1215 
Y0UHES«10.**(XLl/20. ) 

GO  TO  1216 
1215  YOUNESaXLl 
4216  ILHA(J)sYOUNES 
ILMB(J)»-2.*CL1 
FQN(J)=F/1.E9 
UAVL(J)a2.*3.14159/RKl(J) 

TOTRM=R11X  1  ( J )  +RMX2  (  J ) 

T0TRK=RK1(J)+RK2(J) 

TOTRX=XMX 1 ( J ) +XMX 2 ( J ) 

105  CONTINUE 

13  FORMAT  (IX,'  NO  CONVERGENCE  OF  1ST  ITERATION  +  WAVE') 

19  FORMAT  (IX,'  NO  CONVERGENCE  OF  2NO  ITERATION  •••  WAVE') 

139  FORMAT  (IX,'  NO  CONVERGENCE  OF  1ST  ITERATION  -  WAVE') 

199  FORMAT  (IX,'  NO  CONVERGENCE  OF  2ND  ITERATION  -  WAVE') 

DO  4567  IIsl,152 

4567  WRITE(23.3456)FRQ(II),ILMA(II) 

3456  F0RMAT(1X,E12.5,2X,E12.5) 

STOP 

END 

SUBROUTINE  PGRM 

C  RM  SUBPROGRAM  TO  CALCULATE  RADIATION  RESISTANCE 

REALMS  0H,T,F.RKF.RRR,0U,RNMAX,RN.PN,U1,U2,AU1,AU2, 
1RK,RKS:aD,RT,BT,FT,DELT,EXT,A,B,C,D,PC,P1,P3,P21,P22,P2,P,RM1,RM2, 
,RM,X,DET1,DET2,E,G,RTT,RTU,TT,GM,SSECH 
COMMON  F,RKF,RKR,RM1,RM2 
COMMON  O0,T0,B0,RLGH,TU,S0,U0 
COMMON  HI,GA,JCONP,NMAX,FPM 
SSECH(X)»1.-DTANH(X)»DTANH(X) 

DSECH ( X ) =DSQRT ( 1 . -DTANH ( X ) *  DTANH ( X ) ) 

TT=T0+TU 

GM=GA*FPM 

RNsNMAX 

ND3NMAX/2 

RNOsNO 

0H=HI/FPM 

OW»F/GM 

Ul=l.-OH/(OW*OW-OH*OH) 

U2=OW/ ( OW*OW-OH*OH) 

C  Sa+1 

AUlsUl-«-U2-f  1 . 

AU2=Ul-U2-l . 

RKaRKF 

RRS«RK*SO 

RD«RR*DO 

RT=T0»RK 

RTTaTT^RK 

RTU*TU*RK 

BT=0SIN( 0 . 5*RK*B0 ) / ( 0 . 5»RK*B0 ) 

DETlaAUl*TT*DSECH(RTT)»DSECH(RTT)+AU2*DEXP(-2.*RD)*(2.*DO»(Ul+U2-D 
ITANH(RTT) )+TT*DSECH(RTT)»DSECH(RTT) ) 

DET2»AU1 • ( U1-U2+DTANH ( RTT ) ) -AU2*DEXP ( -2 . *RD ) » ( U1 +U2-DTANH ( RTT ) ) 
DBLT»(TO*DSECH(RT)»DSECH(RT)»DTANH(RTU)+TU»DSECH(RTU)»DSECH(RTU)*D 
ITANH(RT) )*DET2+( 1 .+DTANH(RT)*DTANH(RTU) )*DET1 
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EXT=(-1.)*BT/(DELT*RK) 

A=EXT*(AU2*(U1+U2-DTANH(RTU) )*DEXP(-2.*RD)+AU1*(U2-U1-DTANH(RTU) ) ) 

B=EXT*DTANH(RT)*AU1*DSECH(RTU) 

C=EXT*DTANH(RT)*AU2*DEXP(-2.*RD)*DSECH(RTU) 

D»EXT*DTANH ( RT ) »  2 . *U1 *DSECH ( RTU ) 

E=EXT*(AU2*(U1+U2)*DEXP(-2.»RD)+AU1»(U2-U1))*DTAMH(RT)*(-1.) 

G=EXT» ( AUl+DEXP (-2 . »RD ) *AU2 ) *DTAMH( RT) 

PC»3.141592654»F»U0*RK 

Pl=+0 . 5*D*D*DEXP(-2 . *RD) /RK 

P3»0.5*A*A*(DTANH(RT)/RK-T0*DSECH(RT)*DSECH(RT)) 

P21=U1»(B*B*(1.-DEXP(-2.*RD))/RK+C*C*(DEXP(+2.*RD)-1.)/RK-4.*B*C* 

,D0) 

P22=U2* ( B*B* ( 1 . -DEXP ( -2 . *RD ) ) /RK-C*C* ( DEXP ( +2 . *RD )-l . ) /RK) 
P2=0.5*(P21-P22) 

P4=0.5*TU*(E*E-G»G)*SSECH(RTU)+0.5*(E’»E+G*G)*DTAMH(RTU)/RK+E*G*DTA 
1 NH ( RTU ) * DTAMH ( RTU ) / RK 
P=PC*(P1+P2+P3+P4) 

GO  TO  (71,72,93,94,95,96),JCOHF 

71  FT=1. 

GO  TO  99 

72  FT=1. 

DO  920  Nsl.NO 
PNsN 

920  FT=FT+2.*DC0S(PN*RKS) 

GO  TO  99 

93  FT=0. 

DO  930  M=1,ND 
PN»M 

930  FT=FT+2.*DC0S((PN-0.5)*RKS) 

GO  TO  99 

94  FT=0. 

DO  940  N=1,ND 
PN»M 

940  FT=FT+2 . • ( (-1 . ) **N ) •DSIN( ( PN-0 . 5 ) *0 . 5*RKS ) 

GO  TO  99 

95  FT=0. 

DO  950  N»1,ND 
PN=M 

950  FT=FT+2.*DSIM((PN-0.75)*RKS) 

GO  TO  99 

96  FT=0. 

DO  960  Nsl.ND 
PN>N 

960  FT»FT+2.*0SIH((PM-0.25)*RKS) 

99  P»P»FT*FT 
RM1»2.*P 
C  S»-l . 

AU1-U1-U2+1. 

AU2»U1+U2-1 . 

RK>-RKR 

RKS>RK*S0 

RD-RK*DO 

RT=TO»RK 

RTT»TT*RK 

RTU»TU»RK 
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BTaDSIH( 0 . 5*RK*B0 ) / ( 0 . 5*RK*B0 ) 

DETl =AU1 *TT*DSECH ( RTT) *DSECH ( RTT )+AU2*DEXP ( +2 . *RD ) * ( 2 . *D0* < U1 -U2+D 
ITANH(RTT) )+TT*DSECH(RTT)»DSECH(RTT) ) 

DET2=DEXP ( 2 . *RD ) •AU2* ( U1-U2+DTANH ( RTT) )-AUl * ( U1+U2-DTANH ( RTT ) ) 
DELT=(TO*DSECH(RT)*DSECH(RT)*DTANH(RTU)+TU*DSECH(R‘ni)*DSECH(Rnj)*D 
ITANH(RT) )*DET2+( 1 .+DTAHH(RT)*DTANH(RTU) )*DET1 
EXT=BT/(DELT*RK) 

A=EXT*(AU2*(U2-Ul-DTAilH(RTU))*DEXP(+2.*RD)+AUl*(U2+Ul-DTANH(RTU))) 

B=EXT*DTANH(RT)*AU2*DEXP(+2.*RD)»DSECH(RTU) 

C=EXT*DTAMH( RT) *AU1 *DSECH( RTU) 

D=EXT*DTANH  (RT)*2.*'U1*DSECH(RTU) 

E=EXT* ( AU2* ( U2-U1 ) *DEXP ( +2 . *RD ) +AU1 * ( U2+U1 ) ) *  DTAMH ( RT ) • ( - 1 . ) 

G=EXT* ( DEXP ( 2 . 0*RD) *AU2+AU1 ) *DTANH(RT) 

PC=3. 141592654*F*U0*RK*(-1 . ) 

Pls-0 . 5*D*D*DEXP (+2 . *RD) /RK 

P3=0 . 5*A*A» (DTANH(RT) /RK-TO*DSECH(RT)*DSECH(RT) > 
P21=U1*(B*B*(1.-DEXP(-2.*RD))/RK+C*C*(DEXP(+2.*RD)-1.)/RK-4.»B*C* 
,D0) 

P22=U2*<B*B*(1.-DEXP(-2.*RD))/RK-C*C*(DEXP(+2.*RD)-1.)/RK) 

P2=0.5*<P21-P22) 

P4=0.5*TU*(E*E-G»G)*SSECH(RTU)+0.5*(E*E+G*G)*DTAMH(RTU)/RK+E*G*DTA 
1 NH  ( RTU )  *  DTAMH  <  RTU )  /  RK 
P=PC*(P1+P2+P3+P4) 

GO  TO  (81,62.83,84,85,86),JCONF 

81  FT=1. 

GO  TO  89 

82  FT=1. 

DO  820  Nsl.ND 
PN=N 

820  FT=FT+2.*DC0S<PN»RKS) 

GO  TO  89 

83  FT=0. 

DO  830  N»1,ND 
PN=N 

830  FT=FT+2.*DC0S((PH-0.5)*RKS) 

GO  TO  89 

84  FT=0. 

DO  840  N=1,MD 
PN=M 

840  FT=FT+2.*((-l.)»*M)*DSIN((PN-0.5)*0.5*RKS) 

GO  TO  89 

85  FT*0. 

DO  850  Nsl.ND 
PNsN 

850  FTaFT+2.*DSIN((PM-0.75)*RKS) 

GO  TO  89 

86  FT»0. 

DO  860  Nal.ND 
PN=N 

860  FT»FT+2.*DSIM((PH-0.25)*RKS) 

89  P*P*FT»FT 
RI12=2.*P 
RETURN 
END 
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PROGRAM  #2 


C - 

C  PROGRAM  DESCRIPTION: 

C 

C  THIS  PROGRAM  CALCULATES  THE  COMPLEX  DISPERSION  RELATION  FOR 

C 

C  1)  ION  IMPLANTED  MSSV. 

C 

C  2)  UNIMPLANTED  MSSW. 

C 

C  THE  ASSUMED  FORM  IS, 

C 

C  EXP(J(HT-KY)), 

C 

C  CORRESPONDING  TO  A  DISPERSION  FACTOR  OF  THE  FORM, 

C 

C  K  =  BETA  -  J* ALPHA. 

C 

C - 

C 

C  SET  UP  CONSTANTS  AND  DECLARE  COMPLEX  AND  REAL - 

COMPLEX  J,CHO,GAMFO,K3,K4,GK, 

1 GAMMA , E2KU , E2KD , FOFKFO , COSHKV , SINHKV 
REAL  MO 

OPEN(UNITa23 , DEVICE* ' DSK' , DIALOG ) 

PI*3. 1415927 
J=CMPLX(0. 0,1.0) 

GYR0»17.6E+06 

M0*1750.0 

ZEROal.OE-04 

FRACT*1.0E-04 

TYIGIMaO.O 

C 

c - 

C  USER  SUPPLIED  INPUTS 

C - 

TYPE  1200 

1200  F0RMAT( IX, 'ENTER  YIG  AND  SUBST.  THICKNESS  (MICRON)') 

ACCEPT  1201,TYIG,TAL203 

1201  F0RMAT(2F) 

C  YIG  THICKNESS  (MICRONS) - 

C 

C  ALUMINA  SUBSTRATE  THICKNESS  (MICRONS) - 

TYPE  1202 

1202  FORMAT( IX, 'ENTER  H  FIELD') 

ACCEPT  1203, HO 

1203  FORMAT(F) 

C  INTERNAL  BIAS  FIELD  (OERSTEDS) - 

C 

C  GILBERT  LOSS  PARAMETERS  (FO  IN  HERTZ  AND  DELTA-H  MIN  IN  OERSTEDS)- 
F0»3.0E+09 
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DHHINsO.S 

C 

C  DESIRED  FREQUENCY  INCREMENT  IN  OUTPUT  DATA  (MHZ) - 

TYPE  1204 

1204  FORMAT (IX, 'ENTER  FREQ.  INCR.  (MHZ)') 

ACCEPT  1203,DELF 

DELFH*DELF*1.0E6 

C 

C  DESIRED  START  AND  STOP  FREQUENCY  IN  OUTPUT  (GHZ) - 

C  IMPORTANT  ***  :  THE  FSTART  SHOULD  BE  BIGGER  THAN  F(l)  FROM 

C  WU'S  PROGRAM,  AND  FSTOP  SHOULD  BE  SMALLER  THAN  F(152)  FROM  WU'S 
C  PROGRAM. 

TYPE  1205 

1205  FORMAT( IX, 'ENTER  FSTART  FSTOP  (GHZ).  '/ 

NOTE  THAT  FSTART  SHOULD  BE  BIGGER  THAN  FSTART  IN  MSSH' 
PROGRAM  AND  FSTOP  SHOULD  BE  SMALLER  THAN  ITS  FSTOP') 

ACCEPT  1201, FSTART, FSTOP 
FSTSTH=FSTART*1 .0E9 
NUMBER=INT( ( FSTOP-FSTART) » 1000 . /DELF ) 

WRITE(23,111)  FSTSTH,DELFH, NUMBER 
111  F0RMAT(1X,E12.5,2X,E12.5,4X,I5) 

C 

c - ^ - — - 

C 

TYIGIM=TYIGIM*1 .OE-06 
TYIGPU=(TYIG»1 .0E-06)-TYIGIM 
TAL203=TAL203*1 .OE-06 
HSTARTaFSTART*2.0E+09*PI 
WSTOPaFSTOP*  2 . 0E+09*PI 
DELW»DELF»2.0E+06*PI 
WaWSTART-DELtf 
C 

C  MAKE  AN  INITIAL  GUESS  ON  GAMMA - 

GAMF0=CMPLX(0.0,0.0) 

11111  UsU-t-DELU 

F=H/(2.0*PI) 

C 

C  CALCULATE  GILBERT  LOSS  TERMS - - 

DELH=DHMIN» ( FO/F+F/FO) /2 .0 
CHOsCMPLX ( HO , DELH/ 2.0) 

C 

C  CALCULATE  PERMITIVITY  RELATED  CONSTANTS  FOR  FORWARD  PROPAGATION - 

K3=1.0  +  MO/(CHO  -  H/GYRO) 

K4-1.0  +  MO/(CHO  +  W/GYRO) 

C 

C  INCREMENT  NEWTON  RAPHSON  ON  FORWARD  UNIMPLANTED  GAMMA - 

D02  IK2»1, 200,1 

SIMHKV=-J*CSIN(J»GAMF0*TAL203) 

COSHKV»CCOS(J*GAMF0»TAL2O3) 

E2KU«CEXP(2.0*GAMF0»TYIGIM) 

E2KDaCEXP( 2 . 0*GAMF0*TYIGPU) 

FOFKF0= ( K4+1 . 0 ) * ( K3*C0SHKV+SINHKV) *E2KU*E2KD 
1  -(K3-1.0)*(K4*C0SHKV-SINHKV) 

IF( CABS vF0FKF0).LE. ZERO)  GO  TO  1000 
GAMFO«CLOG((K3-1.0)»(K4*COSHKV-SIMHKV) 
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1  /((K4+1.0)*(K3’*COSHKV+SINHKV)))/(2.0*(TyiGPU+TYIGIM)) 

CHREsABS ( REAL ( GAHHA-GAHFO ) ) 

CHIHsABS ( AINAG ( GAHHA-GAMFO ) ) 

OLOREsABS ( REAL ( GAHHA ) ) 

OLDIM-ABS(AinAG(GAHHA) ) 

IF(CHRE.GT.FRACT’»OLDRE)  GO  TO  77 
IF(CHIM.LT.FRACT*OLDIM)  GO  TO  1000 
77  GAMMA=GAMFO 
2  CONTINUE 

WRITE (23, 31)  F.GAMFO.CHRE.CHIH 
31  FORMATdX.'FREQ  =  ME10.3.5X. 'GAMFO  =  '  ,2(E10.3,  IX)  ,4X, 
1/,1X,'CHRE  =  ME10.3,5X. 'CHIM  =  ME10.3) 

STOP 

1000  FREQ=F/1.0E+09 
GK=J*GAMFO 
WRITE(23,101)  GK 
101  F0RMAT(5X,E12.5,2X,E12.5) 

IF(W.LT.USTOP)  GO  TO  11111 

STOP 

END 
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PROGRAM  #3 


THIS  PROGRAM  IS  GENERAL  SYNT.  OF  TRANSDUCERS. 

««««**«*««««*«««***«*«««****««**«*««**««***««[««*«***««** 

**««  IMPORTANT  *in»*»innn»*i»*»*****i»*j»inif***)»)»**** 

A***************************************** ***««*«**«««** 

C  BE  SURE  TO  CHANGE  THE  DIMENSIONS  OF  L.YY.XX  TO  THE 

C  NUMBER  OF  TRANSDUCERS. 

C  TO  SOLVE  THE  MATRIX,  THIS  PROGRAM  MEEDS  A  SUBROUTINE.  IN 
C  THIS  PROGRAM  IT  USES  THE  SUBROUTINE  CALLED  "SIMQ"  WHICH  IS 
C  STORED  IN  THE  MAIN  COMPUETR.  IF  THIS  PROGRAM  IS  USED  ANYWHERE 
C  BUT  "UTA-DEC20  SYSTEM".  YOU  NEED  TO  ADD  MATRIX  SOLVING 
C  PROGRAM. 

DIMENSION  FREQ ( 1 52 } , AAV ( 1 52 ) , FRQ ( I 500 ) , AW ( 1 500 ) . 
1X(1500).L(16),YY(16),XX(16,16),H1(1500),H2(1500) 

DIMENSION  H(1500) 

COMPLEX  K,H,Y1,Y2.Y3,Y4,Y5,X1,X2,X3.X4,X5,Z1,Z2,Z3,K2 
REAL  ISP, LO.LOSSDB.L.PII, KHAR, LOC 
PII=3. 1415927 
TYPE  1000 

C  FOUR  FILE  NAME  IS  NEEDED: 

C  ONE  FOR  READING  F  VS  IL  (LINEAR)  DATA, 

C  ONE  FOR  READING  COMPLEX  K  VALUES. 

C  ONE  TO  STORE  THE  F-IL  DATA, 

C  AND  ONE  TO  STORE  XDUCER  AMD  FILTER  DATA. 

OPEN( UNIT=23 , DEVICE* ' DSK ' , DIALOG ) 

0PEN(UNIT=24 , DEVICE*' DSK' , DIALOG) 

0PEN(UNIT=25, DEVICE*' DSK' .DIALOG) 

0PEM(UNIT*26, DEVICE*' DSK' .DIALOG) 

C  READ  152  DATA  OF  F  VS  I. L( LINEAR) 

READ(23,3)(FREQ(I),AAW(I),I*1,152) 

C  READ  START  FREQ.  ,F  INCREAMENT,  AND  TOTAL  #  OF  DATA. 

REA0(24,333)  FO.DF, NUMBER 
TYPE  1 

C  ACCEPT  #  OF  TRANSDUCERS. 

ACCEPT  2,NUM 
FO«FO-DF 

C  READ  COMPLEX  K  VALUES. 

READ( 24 , 3 ) <K( I ) , 1*1 , NUMBER) 

C  EXPAND  THE  152  F  VS  I.L.  DATA  TO  THE  NUMBER  OF  DATA  FROM 

C  COMPLEX  K  FILE.  THIS  IS  TO  MATCH  F  AND  I.L.  TO  THE  K  VAL 

DO  82  J-l, NUMBER 
FO-FO+DF 
DO  80  1*1,152 
IF(FREQ(I).LE.FO)  GO  TO  81 
DDF1*FREQ< I )-FREQ( I-l ) 

DDAV«AAH(I)-AAH(I-1) 

DDF2*FREQ(I)-F0 

RATI01«DDAH/DDF1 

RESULT*AAW( I )-RATI01 *DDF2 

FRQ(J)*FO 

AH(J)-RESULT 
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60  TO  82 

81  IF(I.EQ.152)  GO  TO  88 

60  TO  80 

88  FRQ(J)=F0 

AV(J)»0.0 
GO  TO  82 
80  CONTINUE 

82  CONTINUE 
TYPE  13 

C  ACCEPT  CEN.  F  IN  HRZ. 

ACCEPT  34,FCENT 

C  FIND  DATA  «  CORRESPONDING  TO  C.F 

DO  90  1=1, NUMBER 
IF(FRQ(I).LE.FCENT)  GO  TO  90 
IFO=I 
GO  TO  91 

90  CONTINUE 

91  CONTINUE 
TYPE  22 

C  ACCEPT  THE  INTEGER  #  (N)  TO  FIND  TRANS.  SPACING  OF  H*LAMBDA 

ACCEPT  34,FSP 

C  CALCULATE  THE  W.L.  OF  C.F. 

XKXK=AIMAG(K(IFO)) 

XXX=2.0»PII*1.0E6/XKXK 

XLAM»INT(XXX) 

S01=XLAM*FSP 

C  CALCULATES  TRANSD.  SPACING  (METERS) 

S0»S01*l.0E-6 
TYPE  1200 

1200  FORMATdX,'  ENTER  I/O  SPACING  (CM)') 

C  ACCEPT  CENT/CENT  SPACING  OF  I/O  TRANSDUCERS  (CM) 

ACCEPT  34, LO 
L0=L0*1.0E-2 

C  FIND  I/O  SPASING  FOR  EACH  TRANSDUCERS. 

DO  6  1=1, NUM 

6  L(I)=L0-KI-1)*S0 

TYPE  5 

C  ACCEPT  B.H  OF  DESIRED  FILTER  (MHZ) 

357  ACCEPT  34, BH 

C  FIND  FREQ.  CORRESP.  TO  THE  B.W  AND  CALCULATING  THE 

C  DATA  ♦  CORRESP.  TO  THESE  FREQS. 

BWW1=(BH*1.0E6)/(2.0*DF) 

BU1=INT(BUV1) 

IF1»IF0~BW1 

IF2aIF0+BWl 

IF(IFl.LT.l)  GO  TO  648 
IF (IF2.GT. NUMBER)  GO  TO  648 
GO  TO  234 
648  TYPE  432 

432  FORMATdX, 'B.W.  IS  TOO  BIG,  TRY  SMALLER  B.W.') 

GO  TO  357 
234  CONTINUE 

RHAR=XKXK*2.0 

C  FIND  W.L.  OF  THE  FIRST  HARMONIC  OF  THE  BAND-PATH  FILTER 

XLAMH-2 . O^PII* 1 . 0E6/KHAR 
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DO  123  I»l, NUMBER 
XKXKsAIHAG(K(I)) 

IF<XKXK.LE.KHAR)  GO  TO  123 
IHAR>1 
GO  TO  321 
123  CONTINUE 

321  IF12=IHAR 

IF(IHAR.LE.IF2)  IF12=IF2+10 

C  EQUATION  FOR  CALCULATING  THE  DESIRED  PHASE (PH=C*F) 

DO  7  Isl, NUMBER 
K2sK(I) 

HHH2s-AIMAG(K2) 

L0C=( ( NUM-1 ) / 2 . ) *SO+LO 

H21=HHH2*LOC 

H22=INT(H21/(2.*PII)) 

H2(I)*H21+H22*2.*PII 
7  CONTINUE 

C  CALCULATE  UPPER  AND  LOVER  VALUES  OF  THE  BAND-PATH  FILTER 

C  IN  THE  LINEAR  SCALE.  (E.G.  IF  HI  IS  1.0  AND  .001,  IT  MEANS 

C  HE  WANT  THE  FILTER  TO  HAVE  O-DB  IN  THE  PASS-BAND,  AND 

C  -60DB  OUTSIDE  THE  PASSBAND. 

DO  20  I»l, NUMBER 
IF(I.GE.IFl)  GO  TO  21 
HKDsO.OOOl 
GO  TO  20 

21  IF(I.LE.IF2)  H1(I)*1.0 

IF(I.GT.IF2)  HKDaO.OOOl 
20  CONTINUE 

DO  25  I:»l, NUMBER 
HHlsHKI) 

HH2=H2(I) 

H(I)aHHl*CMPLX(COS(HH2),SIN(HH2)) 

25  CONTINUE 

TYPE  567 

567  FORMATdX, 'DO  YOU  WANT  SYNT.  ONLY  FOR  THE  FIRST  FUND. 

AHARMONIC  (YES^CR,  N0=1}') 

ACCEPT  2,ICH0 
ISYN«NUHBER 

IF(ICHO.EQ.O)  ISYN«IF12 
DO  30  I1>1,NUH 
Y5»(0. 0,0.0) 

DO  31  Jlal,ISYN 
Y1»K(J1) 

Y2»C0NJG(Y1) 

Y3»Y2*L(I1) 

Y4aCEXP(-Y3) 

31  Y5=H(J1)»Y4»AW(J1)*DF+Y5 

YY(I1)»REAL(Y5) 

30  CONTINUE 

DO  50  Ilal.NUM 
DO  50  I23l,NUM 
X53(0. 0,0.0) 

DO  40  J1«1,ISYN 
Xl-K(Jl) 

X2bC0NJG(X1) 
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X3=X2*L(I1)+K(J1)»L(I2) 

X4«CEXP(-X3) 

40  X5»(AW(Jl)*»2.)*X4*DF+X5 

XX(Il.I2)sREAL(X5) 

50  CONTINUE 

RSsO 

CALL  SIMQ(XX,Yy,NUM,KS) 

Y0=0.0 

DO  100  1=1. NUH 

IF(ABS(YY(I)).GT.YO)  GO  TO  111 
GO  TO  100 

111  YO=ABS(YY(I)) 

100  CONTINUE 

DO  101  1=1. NUH 
YY(I)«YY(I)/YO 

101  CONTINUE 

tfRITE( 26 . 92 ) FRQC 1 ) ,  FRQ<NUHBER) . FCENT. XLAtl, NUH.  SOI . 

ABU . FRQ ( IF 1 ) . FRQ ( IF  2 ) . XLAMH . FRQ ( IHAR ) . FRQ ( 1 ) . FRQ ( I SYN ) . 
B(YY(I).I=1.NUH) 

92  FORMAT( IX, 'NATURAL  CUT-OFF  FRQ.S  ARE:  ',E12.5,'  &  ',E12.5 

A, /'  CENTER  FRQ.  IS:',E12.5,'  WITH  THE  W.L(MIC)  OF  ' 

B, F6.1/,'  #  OF  XDUCERS  =  ',13,  '  WITH  SEPARATION  (MIC)  OF  ' 

C, F6.1/,'  FILTER  B.H  (MH2)=' ,F5.1 ,/'  (START  F=' 

D, E12.5,'  AND  STOP  F=',E12.5,'  )'/ 

E'  H.L  (MIC)  OF  THE  FIRST  HARM*' ,F6. 1 , '  AND  FRQ=' 

F, E12.5,/'  SYNTH.  START  AT  :',E12.5,'  AND  ENDS  AT:  ' 

G, E12.5/'  THE  WEIGHTING  FACTORS  ARE: ' /, 5(1X,E12. 5) ) 

DO  200  I 1=1, NUMBER 
2YY=0 

Z3=(0. 0,0.0) 

DO  201  12=1, NUH 
Z1=-K(I1)*L(I2) 

Z2=CEXP(Z1) 

ZYY=ZYY+(YY(I2)’»*Z.) 

201  Z3=Z2*AH(I1)*YY(I2)+Z3 

Z4=CABS(Z3/ZYY) 

LOSSOB-20 . *AL0G10(Z4 ) 

A1=AIHA6(Z3) 

A2=REAL(Z3) 

A3=A1/A2 

A4«ATAN(A3) 

TBTA»A4*180./PII 
WRITE ( 25 . 3  3 ) FRQ ( 1 1 ) , LOSSDB 
200  CONTINUE 

1  FORMAT( IX, 'ENTER  ♦  OF  TRANSDUERS  ') 

2  FORMAT(I) 

3  FORMAT (2E) 

13  FORMAT( IX, 'ENTER  CENTER  FREQUENCY  (HZ)'  ) 

5  FORMAT (IX, 'ENTER  B.W  (MHZ)  ') 

22  FORMAT( IX, 'ENTER  XDUCER  SPACING  IN  TERMS  OF  LAMBDA') 

26  F0RMAT(1X,5(1X.E11.4)) 

33  F0RMAT(1X,2E) 

34  FORMAT(E) 

333  F0Rt1AT(2E,I) 

1000  FORHATdX, 'THIS  PROGRAM  ASKS  FOR  FOUR  FILE  NAMES.  THE 
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1  FIRST  ONE  IS  FREQ.  VS  I.L.  (152  DATA)  ,  THE  SECOND  FILE 

2  IS  THE  COMPLEX  VALUES  OF  K.,  THIRD  FILE  STORES  THE  F 

3  VS  I.L.{DB),  AND  THE  LAST  FILE  STORES  PARAMETERS.') 

STOP 

END 
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PROGRAM  #4 


C  THIS  PROGRAH  COMBINES  THE  I.L.  DATA  FROM  SINGLE  BAR (ZERO 
C  PATH)  AND  LOOP  (ZERO  PATH).  1/15/85 

OPEM(UMIT=23, devices' DSK' .DIALOG) 

OPEN ( UNIT»24 , DEVICE= ' DSK ' , DIALOG ) 

OPEN (UMIT«25, DEVICE* 'DSK'. DIALOG) 

DO  1  1*1,152 
READ(23.2)F.X 
READ(24,2)F,Y 
RES1*SQRT(X) 

RES2»SQRT(Y) 

RES*RES1*RES2 

1  WRITE(25,3)F,RES 

2  FORHAT(2E) 

3  F0RMAT(1X,2(2X,E12.5))  * 

STOP 

END 
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APPENDIX  3 

"PROGRAM  FOR  CALCULATING  F-I.L  RESPONSE  WITH  EXTERNAL  ATTENUATORS" 


138 


C  THIS  PROGRAM  IS  FOR  HULTIBAR  TRANSDUCER  WITH  EXTERNAL  ATTENUATORES . 
C  OR  THIS  PROGRAM  CAN  BE  USED  TO  CHECK  THE  RESULTS  OF  SYNTH. VALUES 
C  BE  SURE  TO  CHANGE  THE  DIMENSION  OF  L.YY  AND  Y1  TO  THE  #  OF  XDUCER 
DIMENSION  FREQ(152).AAH(152).FRQ(1500).AV(1500). 
1K(1500),L(16),YY(16),Y1(16) 

COMPLEX  K.Z1,Z2.Z3 
REAL  LO.LOSSDB.L.PII 
PII=3. 1415927 
TYPE  1000 

C  PROGRAM  ASKS  FOR  4  FILE  NAME: 

C  FIRST  IS  TO  READ  F  VS  LINEAR  IL  DAT 
C  SECOND  TO  READ  COMPLEX  K  VALUES 
C  THIRD  TO  STORE  THE  F  VS  IL(OB)  OF  FILTER  RESPONSE 
C  AND  THE  LAST  TO  STORE  XDUCER  AND  FILTER  DATA. 
OPEN(UNIT=23.DEVICEs'DSK' .DIALOG) 

OPEN( UNIT=24 , DEVICE* ' DSK ' , DIALOG ) 

OPEN(UNIT=25,DEVICE='DSK'. DIALOG) 

OPEN ( UNIT* 26 . DEVI CE* ' DSK DIALOG ) 
REAO(23.3)(FREQ(I).AAU(I),I*1.152) 

READ(24,333)  FO.DF, NUMBER 
TYPE  1 

ACCEPT  2,NUM 
FO*FO-DF 

READ( 24 , 3 ) ( K( I ) . 1*1 . NUMBER) 

DO  82  J*l. NUMBER 

FO*FO+DF 

DO  80  1*1.152 

IF(FREQ(I).LE.FO)  GO  TO  81 
DDF1*FREQ(I)-FREQ(I-1) 

ODAV*AAU(I)-AAU(I-l) 

DDF2*FREQ(I)-F0 

RATI01*DDAU/D0F1 

RESULT*AAtf ( I )-RATI01 *D0F2 

FRQ(J)*F0 

AV(J)*RESULT 

GO  TO  82 

81  IF(I.EQ.152)  GO  TO  88 

GO  TO  80 

88  FRQ(J)»FO  » 

AV(J)*0.0 
GO  TO  82 
80  CONTINUE 

82  CONTINUE 
TYPE  22 
ACCEPT  34, SOI 
S0*S01*1.0E-6 
TYPE  2000 

2000  FORMAT (IX, 'ENTER  I/O  CENTER/CENTER  SPACING  (CM)') 

ACCEPT  2001,LO 

2001  FORHAT(E) 
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L0«L0»1.0E-2 
DO  6  Isl.NUH 
L(I)»L0+(I-1)»S0 
TlfPE  1200 

1200  FORMAK IX, 'LINEAR  (TYPE  0)  OR  DB  (TYPE  1)  COEF.  ?') 

ACCEPT  1201, NYJA 

1201  FORMATd) 

TYPE  13 

13  FORMAT (IX, 'ENTER  THE  VALUES  OF  THE  ATTENUATORES  ' 

1'  STARTING  WITH  THE  '/,'  TRANSDUCER  CLOSE  TO  INPUT' 

2,'  (F  FORMAT).  ONE  DATA  PER  LINE') 

DO  90  l3l,NUM 
ACCEPT  34,  YKI) 

IF(NYJA.EQ.O)  GO  TO  1203 
YllsABS(Yl(I)) 

Y2»Yll/20. 

Y3»10.0**Y2 

IF(Yl(I).LT.O)  GO  TO  234 

YT(I)=1/Y3 

GO  TO  90 

234  YY(I)=-1./Y3 

GO  TO  90 

1203  YT(I)=Y1(I) 

IF(Y1(I).LT.0)  GO  TO  1205 
Y1(I)=20.*ALOG10(Y1(I)) 

GO  TO  90 

1205  Y1(I)»ABS(Y1(I)) 

Yl(I)a-20.*ALOG10(Yl(I)) 

90  CONTINUE 

34  FORMAT(F) 

HRITE(2.6,92)  FRQ(  1 )  .FRQ(NUMBER)  ,NUM.S01 ,  (Y1  (I) ,  1=1  ,NUM) 
WRITE(26,29)  (YY(I) , 1=1 ,NUM) 

92  FORMAT (IX, 'NATURAL  CUT-OFF  FRQ.S  ARE:  '.E12.5,'  &  ',E12.5/ 

1, '  NUMBER  OF  XDUCERS  =  ',13,'  WITH  SEPARATION  (MIC)  OF  ',F8.2/ 

2. '  VALUES  OF  THE  ATTENUATORES  (IN  DB)  ARE: ' / ,5(1X,E12. 5) ) 

29  FORMATdX,'  OR  THEIR  NORMALIZED  (0  TO  1)  VALUES  ARE:  '/ 

1.5(1X,E12.5)) 

DO  200  I 1=1, NUMBER 

ZYY=0.0 

Z3>(0. 0,0.0) 

DO  201  I2«1,NUM 
Z1»-K(I1)*L(I2) 

Z2«CEXP(Z1) 

ZYY»ZYY+(YT(I2)**2.0) 

201  Z3»Z2*AW(I1)»YY(I2)+Z3 

Z4»CABS(Z3/ZYY) 

L0SSDB-20.*AL0G10(Z4) 

A1«AIMAG(Z3) 

A2«REAL(Z3) 

A3>A1/A2 

A4=ATAN(A3) 

TETA«A4»180./PII 

VRITE(25,33)FRq(Il),L0SS0B 

200  CONTINUE 

1  FORMATdX, 'ENTER  «  OF  TRANSDUERS  ') 
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2  FORHAT(I) 

3  FORHAT(2E) 

22  FORMAT (IX, 'ENTER  XDUCER  SPACING  IN  MICRON  ') 

26  F0RMAT(1X,S(1X,E11.4)} 

33  F0RMAT(1X,2E) 

333  FORMAT(2E,I) 

1000  FORMAT (IX, 'THIS  PROGRAM  ASKS  FOR  FOUR  FILE  NAMES.  THE 

1  FIRST  ONE  IS  FREQ.  VS  I.L.',/,'  (152  DATA)  ,  THE  SECOND  FILE 

2  IS  THE  COMPLEX  VALUES  OF  K.',/,'  THIRD  FILE  IS  FOR  STORING 

3  PROGRAM  INFORMATION  AMD  DATA,  AND  THE  ',/,'  LAST  FILE  IS  TO  STORE 

4  THE  RESULTED  F  VS  I.L.  DATA') 

STOP 

END 
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HSSV  TltAMSVDtSAL  FILTEKS  BASED  ON  CUUtZNT  UEXCBTIHC  IN  NARROW  (10  fia)  TRANSDUCERS 


Y.J.  Acallyan.  J.M.  Omm,  K.W.  R««d.  R.L-  Carctr.  W.A  Davia 


Tha  Cantar  for  Advanead  Elaecron  Davieaa  and  Syataaa  (CARDS) 
Tha  Ualvaralcy  of  Taaaa  at  Arliascen,  Arllttgcon>  Tx.  76019 


ABSTRACT 


wmrm  used. 


A  Nagnetoatatlc  Surface  wave, 
tunable  banoeass  filter  using  currant 
weighted  transducer  arrays  based  on 
transversal  filtering  techniques  has  been 
built.  A  100  hHz  bandeass  filter  was 
realized  with  a  minimum  Insertion 
loss  of  15  dB  and  sidalobe  suedression 
of  ZO  dB.  The  usable  tunabllity 
range  of  the  device  was  from  2  to  3.5  GHz. 


INTROOUCTtON 

Hagnetostatic  Wave  (HSW)  devices 
have  opened  a  new  door  to  high  frequency 
signal  processing  above  I  GHz.  where 
Surface  Acoustic  Wave  (SAW)  devices  have 
high  Insertion  loss  end  are  difficult  to 
fabricate.  One  area  that  has  eluded 
researchers  has  been  tha  develapsNmt  of 
tunable  HSW  transversal  filters.  This 
paper  presents  results  on  the  realization 
of  the  first  synthesized  HSW  transversal 
filters.  Progress  in  this  area  has  bean 
slew,  primarily  as  a  result  of  the 
dispersive  nature  of  HSW' a.  strong 
coupling  of  radiating  eur'rant  elements, 
intaractien  between  neighboring  elements, 
and  reflections  from  the  alaswnts  in  the 
array.  To  overcosw  tha  first  problem,  a 
computer  program  was  developed  to 
synthesize  the  desired  bandpass  filter. 
This  was  based  on  theoretically  calculated 
Insertion  loss  and  freduaney  data  for  a 
zaro-eath  length  delay  line  with  single 
bar  transducers  on  both  the  input  and 
the  output.  The  program  generated 
tha  nonsal  Ized  current  distribution  for 
each  transducer  element.  A  bandpass 
filter  utilizing  nondispersive  waves 
(e.g.  SAW).  ««euld  require  a  currant 
distribution  of  the  form  S1N(X)/X. 
Conversely.  for  HSH's  the  current 
distribution  must  be  distorted  to 
compensate  for  the  dispersion  of  tha  wave. 
In  order  to  reduce  elasMnt  Interaction  and 
reflections  from  the  array,  transducer 
elements  were  made  as  narrow  as  possible 
(10  MS  for  operation  with  a  300  us  canter- 
band  wavelength),  and  In  order  to  limit 
coupling  short  open  circuited  transducers 


THEORY 

Hagnetostatic  waves  (HSW)  are 
basically  magnetically  couoled  dispersive 
waves  that  propagate  I n  a  magnet i ca 1 1 y 
biased  ferrite  material  such  as  Yttrium 
Iron  Garnet  (YIG).  The  theoretical  end 
axparimantal  characteristics  of  the  HSW 
technology  are  well  documented  elsewhere 
Cl. 2. 3. 4]. 

Regardless  of  the  type  of  the  wave 
In  the  propagating  medium  (SAW  or  HSW)  and 
independent  of  the  weighting  technique, 
the  procedure  for  synthes i z i ng  a  f 1 1 ter  i s 
fundamentally  the  same  [5.6].  Por  a 
non-interacting  N  a I ament  transducer  array 
with  a  weighting  factor  of  a  Por  each 
elenmnt  and  path- length  Ig .  the  array 
response  can  be  written  in  tne  formt 


N  -  -K(w)ln  "  * 

6(w)  ■  Z  •nA(u)e  -  Z  fn(u) 

nwl  n*! 


(1) 


where. 

6(e)  ■  array  response 

Me)  ■  e^Ji  (a  and  B  are  loss  factor 
and  t«ave  number,  respectively) 


T  (e)  ■  nth  element  response, 
n 

The  quantity  Me)  denotes  the  complex 
value  of  the  transmission  response  for  a 
slngle-bar/single-bar  transducer  pair  with 
zero  path.  This  transducer  characteristic 
is  Obtained  by  medal  ing  the  transducer 
filaments  as  a  simple  lossless  microstrip 
in  a  manner  dl  scribed  by  Wu  C^l.  A 
pointing  vector  calculation  is  used  to 
define  a  radiation  resistance 
raorasantativa  of  tha  coupling  to  the 
ferrite. 

The  linearized  array  response  given 
I n  equat I on  ( I )  can  be  made  to  fit  the 
^sired  frequency  response  of  the  array. 
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Th«  R.M.s  arror.  e,  can  ba 
ca I cu I area  from  aquat ton  { 2 ) . 


Sa^bing  tna  r^arbtal  dartvabivas  of  bha 
arror  funcblon  wibh  raaaacb  to  aaen  of  bha 
N  paramatars  bo  zaro  yialds  a  mabrtx  of  N 
Hnaar.  Itnaarly  fndaeandanb  aquations, 
for  simollctby,  paramatars  M,  and  Q  ara 
dafinad  ast  'J  < 


1  u,  .  M  -  -K(wjln 
=  “  a„A(ula  } 

ID  I  n»l 


M  -K  (w)l„ 

i(H(u)>£  a^A  (uta  }  du  (2) 

n«l 


(iif  *  II 

Q,  -Raal  /„jA*(w)H(u)a"* 


which  Should  ba  minimizad  with  rasoacb  to 
bha  waighbing  factors  a  : 

m 


Tha  sab  of  aquations  In  (4)  can  ba 

axprassad  as. 


5e 

»4jn  .ii-wi 


'*  N  '*  -K  (ui)t_ 

*(H  anA  (id) a  }dM 

n»l 


w»-*  »K 

t  (u)a  ^ 

,lZ.  .  r"  '  -«*<<«)*„ 

x(R(w)-£  4nA(ula  }du] 


Tha  intagrais  in  aquation  O)  a^a 
conjugata  of  aaen  obhar.  so  bha  sum  can  oa 
axprassao  as  twica  tha  raai  part.  Scsns 
mao I pu I at  I on  y i a I ds • 


"u  ”l2  ”l3  ‘  ■  ‘”lN 

^21  •  . 

*2 

Qz 

X 

s 

a  •  •  •  a  a  a 

• 

a 

Hj  .  ... 

•n 

Qn 

This  sab  of  aquations  can  ba 
algaorabteany  Invartaq  to  obtain  the 
unnorma 1 1 zad  alamant  walgritlng. 

eoaff lei  ants.  onea  the  numarical 
Intagratlon  has  baan  carriad  out  on  tha 
matrix  alamants.  Tha  walghting 

eoaff lei  ants  ware  than  constralnad  to  ba 
real  and  normal  Izaq  bo  unity  <-1  be  'tI). 


Sc  -2  -  -K*(uH_ 

"nr  •  ;::rr'‘^^„,A(MiH(<-)a 


/n/y,|Aluj| 
n*i  • 


-(K  (wltgi  ♦ 

*•  dam  0  (*) 


EAPCRinCNTAt  PROCEOURC 

Tha  16  alamant  filter  in  this 
axaarlmant  was  dasignad  to  have  a  oassband 
of  too  NKz  at  bha  eantar  fraquaney, 
f  "2.7  GHs.  and  as  much  sloaloba 
■ueerasslon  as  bhaorat I ea 1 1 y  pessibla 
using  a  ZS  urn  bhiek  YIQ  film.  Tha 
currant  Is  distrlbutad  among  bha  alamants 
via  a  simpllflad  16  way  pewar-dl vidar. 
which  orovtda  phasa  equalization  for  tarn 
open  cireultad  fingar.  Individual 
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transduesrs  wJthfn  tha  array  ara  all  10  un 
wida  and  4ii¥n  long,  with  only  3rnm  undar 
tha  YIG  (»aa  figura  I  ).  Transducar*  ara 
fad  from  both  ands  with  tha  currant 
waighting  in  aach  transducer  datarminad  by 
tha  placamant  of  an  ooan  cireuitad  gap. 
positionad  ralativa  to  tha  longitudinal 
array  axis.  As  an  axampla.  a  gap 
at  tha  cantor  of  a  transducar  would 
rasult  In  a  zaro  nat  currant 
contribution  for  that  alamant.  It  was 
assumad  that  in  tha  caso  of  an 
opon-c i rcu I ted  mlcrostrlp  of  length 
l  <<  x/4.  with  fringe  affects 
neglected,  the  current  exhibits  linear 
variation  along  the  length.  Therefore, 
the  gap  in  each  alamant  was 
positionad  so  that  the  average  current  on 
opposing  elements  had  tha  raouirad 
sign  and  value  dictated  by  the  synthesis 
program.  In  order  to  reduoa  the  harmonic 


rosponses  of 
spaced  loop 


the  array, 
transducer 


Is 


short ad 
usad  on 


i/2 

the 


output. 

Tha  transducers  ara  fabricated  from 
up-piatad  gold  on  10  mil  thick. 
I"xl"  alumina  substratas.  The  3nfn  wide. 
25  iMi  thick  YIG  film  is  angla-lapoad  at 
tha  ands  to  minimize  reflections  there. 


EXPERIMENTAL  RESULTS 

Figure  2  shows  tha  theoretical 
pradiction  of  tha  insertion  loss  versus 
fraouancy  for  a  typical  filter  with  tha 
desired  characteristics  described  earlier. 
Tha  minimum  Insertion  loss  of  the 
fabricated  device  was  observed  to  be 
-2SdB.  which  is  mainly  due  to  the 
power-divider  mismatch  and  associated  20:1 
VSWR  i n  the  passband  of  the  f 1 1 tar . 
Figure  3  is  the  transmission  response  for 
tha  dev ice  In  the  frequency  range  of  1.5 
GHz  to  3.0  GHz.  In  this  picture  the 
output  of  the  device  Is  amplified  to  show 
the  side  lobe  suppression,  otherwise  not 
clearly  distinguishable  from  direct 
breakthrough . 


Os 


Figure  2:  Tbeertclcsl  prtdletlen  of  S^2* 


Figuru  1:  Clreulc  pscesm  el  the  100  MHs  cuseble 
baadpess  fllcer.  The  pesicies  ot  the 
gap  is  aarkad  with  at  cha  laic  aida 
ler  aach  craaadnear. 

Actual  dlaanaion:  1"X1". 
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Th«  tnaorecical  simulation  the 
tunabIHty  of  tha  filter  Is  shown  In 
f < gura  4 .  This  was  dona  by  sucar I mpos I ng 
tha  calculated  transmission  rasoonses  of 
tha  filter  at  four  different  bias 
fields.  Figure  5  shows  tha  experimental 
results  obtained  by  varying  tha 
magnetic  bias  field.  Figures  4  and  -5 
clearly  demonstrate  a  decrease  in  the 
oaasband  for  higher  center  f reouenc I  as . 
This  is  due  to  the  band  limiting 
characteristics  of  MSSW  delay  line 
dictated  by  the  applied  bias  field. 


Figurs  4;  Thaarseleal.  vecClctlen  tot  twubUity 
ot  Sj2. 


2.0  O.IS  CHs/div  3.S 


CHs 

Figure  3:  Tunsblllcy  of  S^2  (200  -  430  Os). 


CONCLUSION 

The  current  weighting  technique 
using  open  gap  transducers  to  realize  HSSU 
filters,  demonstrates  a  good  corre I at I on 
between  tha  theoretical  predictions  and 
obtained  performance.  With  this 
taehn I qua •  it  Is  poss I b I e  to  construct  a 
tunable  bandpass  filter  with  sidalobe 
suppression  as  high  as  20  dB.  with  less 
than  3  dB  amp 1 1 tuda  ripple  .  .  Tha 
vriation  of  the  minimum  I nsart i on  I oss  in 
the  tunable  range  of  tha  device  can  be  as 
low  as  2  dB.  An  improvement  of  about  10 
dB  on  the  minimum  insertion  loss  was 
achieved  by  a  narrow  band  matching  of  tha 
existing  power-divider.  Therefore,  we 
be I i eve  further  I mprovament  in  deli ver I ng 
the  maximum  power  to  each  finger  via 
construction  of  a  better  matched  broad 
band  power-divider  can  bring  the  minimum 
insertion  loss  close  to  -10  dB. 
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Technical  Feature 


Delays  using  Magnetostatic 
wave  Technology 
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Variable  time  delays  are  necessary  in  phased  array  systems  to  prevent  phase  spuinting  and  pulse  stretching. 
Methods  tor  providing  these  time  delays  include  an  assortment  of  fixed  cables,  ferrite  loaded  cables,  surface  acoustic 
wave  iSAW)  devices  and  magnetostatic  wave  iMSW)  devices.  Fixed  cables  are  bulky,  limiting  the  number  that  can 
be  employed  per  system.  Fernte  loaded  cables  and  SAW  devices  are  applicable  primanly  at  frequencies  below 
1  GHz  and  provide  relatively  small  delay  differentials.  MSW  wave  technology  is  capable  of  operating  at  frequencies 
up  to  20  GHz  and  providing  differential  time  delays  on  the  order  0/  tens  of  nanoseconds.  An  MSW  device  has  recently 
been  demonstrated  with  a  bandwidth  greater  than  200  MHz  centered  at  3  GHz.  This  device  has  a  phase  error 
across  the  band  as  low  as  ^  and  is  capable  of  providing  nearly  50  nS  differential  delay.  Thus.  MSW  technology 
appears  to  be  the  most  promising  technique  for  the  next  generation  of  phased  array  systems. 


Introduction 

The  need  for  enhanced  perform¬ 
ance  and  miniaturization  in  ohased 
array  system  components  has 
fueled  the  exploration  of  new  tech¬ 
niques  for  providing  time  delays. 
Ideally,  a  time  delay  component 
snouid  tie  small  ana  rapidly  tunable 
over  a  delay  range  of  a  few  tens  of 
nanoseconds,  should  have  excel¬ 
lent  phase  linearity  characterisitics. 
and  Should  be  inexpensive,  in  this 
paper  we  describe  recent  advan¬ 
ces  in  the  state  of  the  art  for  one  new 
tunable  time  delay  technique,  the 
magnetostatic  wave  iMSViO  delay 
line. 

True  time  delays  perform  two 
principal  functions  m  phased  array 
systems:  they  eliminate  "phase 
squinting '  m  broadband  beams  and 
they  allow  the  undistorted  trans¬ 
mission  and  reception  of  narrow 
pulses.  Consider  the  first  problem. 
In  a  phased  array  antenna  the  beam 
IS  steered  by  adjusting  the  phase  of 
the  electromagnetic  signals  from 
each  of  a  large  number  of  radiating 
elements  so  that  the  radiated  waves 
add  coherently  only  m  a  specified 
direction  if  frequency  independent 
phase  shifters  such  as  diode 
phasers  are  used  to  provide  the 
needed  phase  shift  then  for  a  given 
steering  angle  the  signal  will  be 


strictly  coherent  only  at  .one  fre¬ 
quency. 

Consider  an  array  of  radiators 
with  spacing  d  as  shown  m  Figure 
1  If  the  steering  angles  is  i  and  0  is 
the  path  difference  between  adja¬ 
cent  raoiators.  then  the  relation 
oetween  n  aro  D  is  given  by 
d''sin(»)-*0.  Thus,  tr’e  'eouired 
phase  at  each  'aciator  for  wave- 
length  L  will  be  2’-’  0  L  C  eariy. 
a  frequency  moepencent  pnase 
shifter  can  satisfy  this  reouirement 
at  only  one  frequency.  Other  fre¬ 
quency  components  will  sum  co¬ 
herently  in  slightly  different  direc¬ 
tions.  and  the  net  system  result  will 
be  high  sidelobesoutsioe  a  relative¬ 
ly  narrow  passband.  This  is  illustrat¬ 
ed  in  Figure  2  by  a  simuiatea  beam 
pattern  from  an  array  of  radiators. 

For  this  calculation,  a  uniform  lin¬ 
ear  array  with  10  isotropic  elements 
was  assumed  (i.e.  no  attempt  was 
made  to  modify  the  sintxi/x  re¬ 
sponse  of  the  passband).  The  de¬ 
sired  beam  angle  for  the  array  was 
4S'.  the  element  spacing  was  S 
wavelength  and  the  phase  error 
was  simulated  using  a  random 
number  generator.  Calculations 
were  made  for  two  center  frequen¬ 
cies  separated  by  250  MHz  (2.875 
and  3.125  GHz),  with  the  assump¬ 
tion  that  beamstcering  was  accom¬ 


plished  using  ideal  phase  shifters. 
Figure  2  shows  clearly  that  the  com¬ 
mon  passband  of  these  two  fre¬ 
quencies  has  been  narrowed  by  the 


-  •  aeann  it»ere<J':rro'.,gn  .irqie  0 


beamsteering.  In  effect,  a  broad¬ 
band  imaging  system  has  been 
squinted"’  mto  a  narrower-bano 
one. 

The  second  problem,  the  pulse 
distortion  of  a  narrow  pulse,  is  inde¬ 
pendent  of  bandwidth  and  will  anse 
whenever  the  steering  angle  is 
large.  This  is  illustrated  m  Figure  3. 
Let  Ta  represent  the  time  required 
for  an  electromagnetic  signal  to 
travel  across  an  antenna  array.  If  the 
beam  is  steered  an  angle  i  away 
from  the  array  normal,  then  the  lead- 
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Fig.  2  Beam  snap*  for  an  ideal  pnasa  shiftar  at  two  frequancias 
saparatod  by  2S0  MHz. 


Fig.  3  Pulsa  distortion  dua  to  uncompansated  tima  oaiay  across 
tna  array. 


ing  edge  of  a  transmitted  pulse  from 
the  near  side  of  the  array  will  arrive 
at  a  target  TA*sin(^  earlier  than  the 
leading  edge  of  the  pulse  trans* 
mitted  from  the  far  side  of  the  array. 
In  effect  the  pulse  width  will  be 
stretched  by  an  amount  TA*sin(a). 
Since  Ta  is  usually  on  the  order  of 
a  few  tens  of  nanoseconds,  the  ef* 
feet  is  important  only  for  narrow 
pulses.  However,  narrow  pulses  are 
mandatory  when  high  resolution 
imaging  of  a  target  is  required,  and 
thus  the  problem  is  highly  signifi* 
cant 

Fixed  Delay  Eements 
In  most  present  day  systems,  time 
delays  are  provided  by  an  assort¬ 
ment  of  fixed  transmission  lines,  ei¬ 
ther  stnpline  or  cable.  An  excellent 
elementary  review  of  this  approach 
has  been  given  by  BrooKner.'  For 
large  angle  sweeping,  a  3-bit  phase 
shifter,  consisting  of  three  unequal 
len^  stnplines.  is  provided  at  each 
radiating  element  These  lines  have 
paths  equivalent  to  1/2.  1  /4.  and 
1  /8  wavelength  and  can  be  com¬ 
bined  to  give  phase  differences  of  0* 
to  360*  in  steps  of  45*.  The  three 
stnplines  together  are  less  than  one 
waveien^  long,  so  that  the  total 
phase  shift  is  never  more  than  360*. 
Thus,  although  broadband  opera¬ 
tion  IS  realized,  high  resolution,  nar¬ 
row  pulses  a  few  tens  of  wave¬ 
lengths  wide  cannot  be  recon¬ 
structed  unambiguously  using 
these  elements  alone. 


For  narrow  pulse  work,  the  array 
is  commonly  divided  into  subarrays, 
and  each  of  these  units  is  provided 
with  a  3-bit  phase  shifter  comprised 
of  tong  cables.  For  a  system  such  as 
the  COBRA  DANE  radar,  delays  as 
long  as  64  wavelengths  can  be  ob¬ 
tained  for  each  subarray.  Although 
this  technique  is  believed  to  be  ade¬ 
quate  for  most  current  applications, 
there  is  necessarily  some  distortion 
of  the  pulse  due  to  the  finite  delay 
across  each  subarray.  However,  it  is 
clearly  impractical  to  provide  each 
element  of  the  array  with  long  delay 
lines,  since  something  like  1 00  miles 
of  cable  would  be  required  for  a 
radar  like  COBRA  DANE  Even  us¬ 
ing  the  subarray  approach,  about 
one  mile  of  cable  is  necessary.  It 
would  clearly  be  desirable  to  re¬ 
place  these  great  lengths  of  cable 
with  compact  modules:  therefore, 
considerable  effort  has  been  ex¬ 
pended  in  recent  years  in  explonng 
methods  of  varying  delays  electron¬ 
ically. 

The  principal  approaches  inves¬ 
tigated  to  date  include  ferrite  loaded 
helical  transmission  iines.^  surface 
acoustic  wave  (SAW)  delay  lines^'^ 
and  magnetostatic  wave  (MSW)  de¬ 
vices.  The  delay  is  changed  in  the 
femte  loaded  line  by  applying  an 
external  magnetic  field.  This 
changes  the  permeability  of  the  fer- 
nte.  the  impedance  of  the  trans¬ 
mission  line  and.  hence,  the  velocity 
of  the  electromagnetic  wave.  The 
velocity  of  the  surface  acoustic 


wave  is  changed  by  applying  a  large 
voltage  across  the  delay  path.  Both 
the  fernte  loaded  line  and  the  SAW 
delay  line  operate  at  low  frequen¬ 
cies  (below  2  GHz)  and  provide  rel¬ 
atively  small  differential  delays  (<20 
nS).  In  contrast  magnetostatic  wave 
(MSW)  devices  have  the  potential  of 
operating  at  frequencies  up  to  20 
GHz  with  differential  delays  greater 
than  50  nS.  Thus,  this  paper  is  con¬ 
cerned  with  the  current  status  of 
MSW  technology. 

Magnetostatic  Delay  Line 

Considerable  progress  nas  been 
made  recently  toward  the  realiza¬ 
tion  of  an  MSW  variable  time  delay 
device  that  can  operate  in  the  mi¬ 
crowave  frequency  range  and  pro¬ 
vide  electronically  tunable  delays  of 
a  few  tens  of  nanoseconds.  This 
device  has  been  developed  under 
the  sponsorship  of  Rome  Air  Devel¬ 
opment  Center  through  a  loint  effort 
of  Rockwell  International.  The  Uni¬ 
versity  of  Texas  at  Arlington  and 
North  Carolina  State  University. 

The  propagating  medium  for 
magnetostatic  wave  iMSW)  devices 
consists  of  an  epitaxial  fernte  film 
(usually  yttrium  iron  garnet  or  YIG) 
on  a  garnet  substrate.  This  structure 
is  usually  prepared  m  the  form  of  a 
bar  with  input  and  output  trans¬ 
ducers  at  each  end.  The  device  is 
placed  between  the  poles  of  a 
magnet  and  an  RF  signal  is  fed  into 
one  transducer.  An  illustration  of  a 

[Connnutd  on  pago  1 1 2| 
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pig.  4  Excitanon  ot  MSSW.  MSBVW  and  MSP/W  on  a  YiG  daiay  ime  using  a  snorted 
4i«ctroce  transducer,  [from  Sbgiiu  and  Setnares°1 


Standard  MSW  delay  une  is  snown 
in  Figure  4.®  At  a  specific  comoma- 
tion  of  signal  frequency  and  mag¬ 
netic  field  strength,  magnetic  som 
waves  will  be  launched  from  the 
input  transducer  down  the  bar. 
These  waves  are  recon veneo  into 
electromagnetic  energy  at  the  out¬ 
put  transducer.  Since  the  velocity  of 
magnetostatic  waves  m  the  YIG  is 
some  three  orders  of  magnitude 
smaller  than  electromagnetic 
waves  in  free  space,  a  substantial 
delay  is  realized  with  a  delay  path 
of  about  1  cm.  Magnetostatic  waves 
are  somewhat  analogous  to  surface 
acoustic  waves,  but  they  have  two 
advantages  for  phased  array  sys¬ 
tems  applications;  they  operate  from 
1  to  20  GHz.  and  their  frequency  of 
operation  can  be  electronically 
tuned  by  changing  the  value  ot  the 
magnetic  bias  field. 

MSW  delay  lines  are  inherently 
dispersive,  i.e.  the  group  velocity 
does  not  in  general  equal  the  phase 
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/eiocity.  The  scecific  dispersive 
cnaractenstics  of  an  MSW  delay 
re  are  oecenoent  on  a  number  ot 
/anabies.  but  one  useful  feature  is 
that  the  siooe  of  the  dispersion  can 
be  either  oositive  or  negative,  de- 
oending  on  the  orientation  of  the 
magnetic  bias  field.  This  aspect  nas 
been  applied  m  the  oesign  of  a  var- 
abie  delay  .'me  that  operates  at  the 
center  frequency  near  3  GHz  and 
exhibits  a  oeiay  differential  of  more 
man  40  nS. 

A  schematic  diagram  of  the  ap- 
oroacn  is  shown  m  Figure  5.  Two 
MSW  delay  imes  navmg  dispersive 
characteristics  with  slopes  of  oppo¬ 
site  Signs  are  cascaded  together,  if 
the  oispersions  are  linear  and  the 
slopes  are  equal  m  absolute  mag¬ 
nitude.  then  the  net  dispersion  will 
be  zero  across  the  passband.  if  one 
delay  line  is  provided  with  a  con¬ 
stant  bias  while  the  bias  on  the  other 
IS  varied,  the  net  delay  will  change 
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soacirg  between  the  film  and  the 
ground  plane.  For  the  M.SSW  delay 
line  a  variable  grcura  oiane  «as 
prcv^deo  Ti-ese  technaues  are  de¬ 
scribed  in  mere  detail  oy  Chang  et 
al*  and  In  rere'-erce  3 
Tre  ddc^agTg  deserves  some 
dommert.  it  :s  moortant  that  the 
f’eids  seen  by  the  delay  imes  be 
uniform  to  orevent  additional  un- 
-varted  frecuency  disoersion  Such 
uniformity  s  readily  achieved  with 
arge  aboratory  magnets,  but  small 
oacKagea  components  -suaily  pro¬ 
duce  non-uniform  magnetic  fields. 
Our  approach  was  to  place  the  bi¬ 
asing  magnets  at  opposite  ends  of 
the  oackages  and  focus  their  fields 
along  the  central  axis  by  means  of 
a  senes  of  smaller  magnets  with 
opposing  fields  suitably  spaced 
along  the  remaining  sides,  top  and 
bottom.  As  the  fields  required  were 
relatively  small  t600  oe).  ferrite  mag¬ 
nets  were  used  throughout  In  the 
photograph  the  magnets  are  the 
dark  colored  rectangles,  while  the 
light  rectangles  are  aluminum  spac- 
[Continuad  on  pag»  1 16| 
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but  remain  flat  as  a  function  of  fre¬ 
quency  over  the  passbano  yVe  will 
-■efer  to  a  device  of  tn.s  tyce  23  a 
MSW  cascaded  delay  me  C3Li  a 
laboratory  version  of  the  aoorcacn 
was  first  reported  by  Sethares  2ra 
Owens,  and  an  improveo  acora- 
tory  version  was  reported  by  me 
authors  of  this  article  m  1984  *  This 
device  has  now  been  packaged  and 
is  sfiOwn  in  Figure  6. 

The  negatively  sloped  dispersion 
IS  provided  by  a  backward  volume 
wave  (MSBVW)  delay  ime  while  a 
surface  wave  iMSSW)  delay  une 
gives  the  positive  slope  character¬ 
istic.  The  MSBVW  mode  is  obtained 
when  the  bias  field  is  parallel  to  the 
direction  of  propagation  The 
MSSW  mode  is  produced  when  the 
field  IS  applied  m  the  plane  of  the  fifm 
and  IS  perpendicular  to  the  direction 
of  propagation.  One  problem  inher¬ 
ent  in  MSW  dispersion  is  that  it  is  m 
general  non-linear.  Thus,  some 
method  for  linearizing  these  char¬ 
acteristics  must  be  used.  For  the 
MSBVW  defay  line,  the  dispersion 
was  linearized  by  adjusting  the 

149 


(From  page  114|  AOKINS 


Soiitron/MjGrowave 


’suracniE* 

maHTWBU 


SF2912-6605 


Piug/JacK  Adaoter 
SF2994^2 


Fiaige  Mount  • 
Plug  Peceotacte 
SF2g61-6601 


SoUtron/Micrawaves  new  l8GHz 
SMA  Cubed  Right  Angle  Connectors 
are  designed  with  a  unioue  cubed 
body,  teflon  insulators  and  a  one 
piece  pre-bent  center  contact.  The 
result?  Lower  VSWR  and  insenion 
loss  at  higher  treouencies!  Also, 
these  high  frequency  Super  Cube 
Connectors  are  available  at  a  much 
lower  cost  than  standard  Radius 
Right  Angle  types. 

Typical  specifications  include: 

•  FREQUENCY  RANGE; 

DC  to  I8GH2 

•  '/SWR;  1.05  +  008  x  fGH? 

•  INSERTION  LOSS; 
OSdBxVTSHl 

Our  cubed  Right  Angle  Cormectors 
are  available  m  many  configurations 
including  receptacle,  cabled  and 
adapter  types. 

Contact  us  today  tor 
complete  intotmatian  and  literature. 


(Solitron/Microwave] 


►  '  M 


Fig.  6  Packaged  CDL 


Fig.  T  Linear  disperaions  with  MSSWV  ara  MSSW  ceiay  imes. 
lai  MS8VW;  bl  MSSW 


ers.  The  uniformity  of  the  biasing 
fields  was  excellent  using  this  tech¬ 
nique.  The  variable  field  was  pro¬ 
vided  by  means  of  coils  in  the 
MSBVW  package,  and  the  outer 
shell  of  both  components  was  stain  - 
less  steel. 

The  frequency  dispersive  char¬ 
acteristics  of  the  packaged  MSSW 
and  MSBVW  delay  lines  are  shown 
in  Figure  7.  and  the  net  frequency 
dispersion  from  the  complete  COL  is 
given  in  Figure  8.  The  figure  shows 
a  deiay  differentiai  of  34  nS  (approx¬ 
imately  too  wavelengths  at  center 
frequency)  when  the  current 
through  the  coil  is  changed  from  0 
to  .6  ampa  The  maximum  delay 
change  that  has  been  demonstrated 
with  this  configuration  is  47  nS. 


From  a  visual  inspection  of  these 
^lustrations,  the  linearity  of  the  indi¬ 
vidual  delay  lines  and  the  flatness  of 
the  cascaded  device  appear  to  be 
excellent.  For  a  more  quantitative 
measure  of  these  characteristics, 
the  root  mean  square  (RMS)  of  the 
phase  oeviation  from  linearity 
across  the  band  was  computed. 
This  calculation  was  based  on  the 
fact  that  the  phase  is  linear  across 
a  non -dispersive  passband.  Phase 
data  as  a  function  of  frequency  were 
obtained  from  the  automatic  net¬ 
work  analyzer  and  were  fitted  to  a 
linear  curve.  The  deviation  of  the 
measured  phase  from  the  calculat¬ 
ed  curve  was  taken  to  be  the  phase 
error,  and  the  RMS  value  was  com¬ 
puted  from  these  data  using  stan- 
[Continucd  on  pago  118| 
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Fig.  8  Oisporsion  curvM  tor  tn«  packagod  COL  at  mraa  currant  settings. 


dard  statistical  techniques.  The  crit¬ 
ical  charactenstics  of  the  device  are 
given  in  the  Table  1. 

The  characteristics  listed  above 
come  very  close  to  meeting  system 
needs.  Refernng  back  to  the  discus¬ 
sion  on  fixed  delay  lines.  COLs  with 
these  characteristics  could  replace 
the  assortment  of  tong  cables  at  the 
subarray,  reducing  a  mile  of  cable 
to  a  set  of  compact  (<1 20  cm^  com¬ 
ponents.  However,  it  still  will  be  nec¬ 
essary  to  provide  the  radiators  with 
3-bit  0*  to  360*  phase  shifters.  Since 
the  MSW  delay  lines  are  frequency 
dispersive,  the  group  velocity  line¬ 
arized  by  the  COL  is  not  equal  to  the 
phase  velocity.  Thus,  the  COL  can¬ 
not  simultaneously  provide  a  given 
absolute  phase  and  the  required 
group  delay.  To  provide  phase  co¬ 
herence  across  the  array.  0*  to  360* 


phase  shifters  will  be  necessary. 
This  should  not  be  a  serious  draw¬ 
back.  however,  since  these  phase 
shifters  are  reasonably  small  and 
can  be  integrated  on  a  semi¬ 
conductor  chip. 

Simulated  beam  shapes  using 
COLs  plus  phase  shifters  are  shown 
in  Figures  3  and  tO.  Except  for  the 
addition  of  MSW  delay  lines  to  the 
beamsteering  mechanism,  the  as¬ 
sumptions  here  are  identical  to 
those  used  in  calculating  Figure  2. 
As  before,  the  beam  shape  is  calcu¬ 
lated  for  two  center  frequencies: 
2.875  and  3.125  GHz.  For  the  two 
figures  here,  the  delay  lines  have 
RMS  phase  errors  of  8*  and  12*. 
repsectively.  It  is  clear  that  beam 
squinting  has  been  dramatically  re¬ 
duced  using  components  with 
these  charac^stlcs. 
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Conclusion 

As  a  variable  time  delay  compo¬ 
nent  the  MSW  COL  has  two  clear 
advantages  over  competing  tech¬ 
nologies  such  as  ferrite  loaded  ca¬ 
bles  and  surface  acoustic  wave  de¬ 
lay  lines:  the  frequency  of  operation 
IS  well  into  the  microwave  region  (3 
GHz  for  the  present  device),  and 
delay  differentials  on  the  order  of 
tens  of  nanoseconds  (47  nS  for  the 
present  device)  are  obtainable.  With 
respect  to  the  first  feature,  note  that 
although  the  present  device  has  a 
center  frequency  of  3  GHz.  MSW 
technology  is  capable  of  performing 
at  frequencies  up  to  20  GHz.  This  is 
in  contrast  to  surface  acoustic  wave 
devices  that  are  essentially  con¬ 
fined  to  a  frequency  regime  of  1  GHz 
or  less  and  the  ferrite  loaded  helix 
that  operates  at.  or  below  the  VHF 
region  (300  MHz).  The  differential 
delay  obtainable  with  the  MSW  de¬ 
vice  also  is  impressive.  A  delay 
change  of  47  nS  has  been  obtained 
with  the  MSW  COL  while  a  maxi¬ 
mum  variation  of  4  nS  has  been 
reported  for  the  helix  and  about  10 
nS  for  the  SAW  delay  line.  In  the 
SAW  case,  a  voltage  on  the  order  of 
kilovolts  was  required  to  realize  the 
reported  change. 

The  main  drawbacks  to  the  MSW 
device  are  sensitivity  to  temperature 
Change,  moderate  insertion  loss 
and  low  oower  handling  capability 
The  power  handling  capability  is  in¬ 
trinsic  to  the  magnetostatic  mterac- 
t  on  process,  but  me  insertion  loss 
c.'iaracteristic  and  temperature  sta¬ 
bility  can  be  improved  with  further 
research  Thus,  the  results  reported 
here  indicate  that  MSW  devices  will 
be  capable  of  meeting  many  of  the 
time  delay  requirements  of  future 
generation  pnased  array  systems. 
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